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This article presents detailed discussions and calculations of the recent paper "Quantum 
Regge calculus of Einstein-Cartan theory" in Phys. Lett. B682 (2009) 300. The Euclidean 
space-time is discretized by a four-dimensional simplicial complex. We adopt basic tetrad 
, and spin-connection fields to describe the simplicial complex. By introducing diffeomorphism 

, and local Lorentz invariant holonomy fields, we construct a regularized Einstein-Cartan 

theory for studying the quantum dynamics of the simplicial complex and fermion fields. 
This regularized Einstein-Cartan action is shown to properly approach to its continuum 
counterpart in the continuum limit. Based on the local Lorentz invariance, we derive the 
_C dynamical equations satisfied by invariant holonomy fields. In the mean-field approximation, 

we show that the averaged size of 4-simplex, the element of the simplicial complex, is larger 
^ ' than the Planck length. This formulation provides a theoretical framework for analytical 

£f} calculations and numerical simulations to study the quantum Einstein-Cartan theory. 
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Since the Regge calculus 

HQ 

was proposed for the discretization of gravity theory in 1961, 
many progresses have been made in the approach of Quantum Regge calculus [j, |J] and its vari- 
ant dynamical triangulations [5]. In particular, the renormalization-group treatment is applied 
to discuss any possible scale dependence of gravity [3|. Inspired by the success of lattice regu- 
larization of non-Abelian gauge theories, the gauge-theoretic formulation [0\ of quantum gravity 
using connection variables on a flat hypercubic lattice of the space-time was studied in the La- 
grangian formalism. The canonical quantization approaches to the Regge calculus in Hamiltonian 
formulation are studied in Ref. {?]]. A locally finite model for gravity has been recently proposed 
8(]. All these studies are very important steps to understand the Einstein general relativity for 
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based on the 



gravitational fields in the framework of quantum field theory. In the brief paper 
scenario of quantum Regge calculus, we present a diffeomorphism and local Lorentz invariant (i.e., 
local gauge-invariant) regularization and quantization of Euclidean Einstein-Cartan (EC) theory. 
Detailed calculations and discussions are presented in this article. 

The four-dimensional Euclidean space-time is discretized by a simplicial complex, analogously 
to the formulation of the Regge calculus. In the framework of the Einstein-Cartan theory, we 
adopt basic gravitational variables, i.e., a pair of tetrad and spin-connection fields to describe 
the simplicial complex. Introducing diffeomorphism and local Lorentz invariant (i.e., local gauge- 
invariant) holonomy fields in terms of tetrad and spin-connection fields along loops, we propose 
an invariantly regularized EC theory for the dynamics of simplicial complex, which couples to 
fermion spinor fields. We show that in the continuum limit when the wavelengths of tetrad and 
spin-connection fields are much larger than the Planck length, this regularized EC action properly 
approaches to the continuum EC action. The quantum dynamics of the simplicial complex is 
described by the Euclidean partition function that is a Feynman path-integral overall quantum 
tetrad, spin-connection, and fermion fields with the weight of regularized EC action. Based on 
local gauge invariance, we derive the dynamical equations satisfied by invariant holonomy fields of 
tetrad, spin-connection, and fermion fields. In the mean-field approximation, we show the averaged 
size of 4-simplex (and its 3-simplex and 2-simplex), elements of the simplicial complex, has to be 
larger than the Planck length. This formulation provides a theoretical framework for analytical 
calculations, in particular, numerical simulations to study the Einstein-Cartan theory as a quantum 
field theory. 

This article is organized as follows: In Sec. HH we give a brief review of the continuum EC theory. 
In Sec. IIIIl we discuss the regularized EC theory based on (1) the description of simplicial complex 
by tetrad and spin-connection fields; (2) parallel transport equations in simplicial complex; (3) 
invariant holonomy fields and regularized EC action and their continuum limit; (4) the Euclidean 
partition function. In Sees. IIVI and W \ we study chiral gauge symmetric bilinear and quadralinear- 
fermion actions, and derive dynamical equations for holonomy fields. In Sec. IVH we adopt the 
method of the mean-field approximation to show the averaged size of the 4-simplex has to be 
larger than the Planck length. In the last section, we give some concluding remarks, and detailed 
calculations are arranged in Appendices [A] |Bl [Cj [D] and [El 
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II. CONTINUUM EINSTEIN-CARTAN THEORY 

The basic gravitational variables in the Einstein-Cartan theory constitute a pair of tetrad and 
spin-connection fields [e M a (x), lo'^(x)], whose Dirac-matrix values 

e p (x) = e A1 ct (x)7 a and = u ab (x)a ab . (f) 

The fields e p {x) and are 1-form real fields on the four-dimensional Euclidean space-time 

1Z , taking values, respectively, in the local Lorentz vector space Vc and in the Lie algebra so(4) 
of the Lorentz group SO(A) of the linear transformations of Vc preserving 5 ab = (+,+,+,+). In 
this local Lorentz vector space Vc, fermions are spinor fields ip(x), Dirac 7 matrices obey 

{la,lb} = ~25 ab , (2) 

7T = — 7 a and 7^ = — 1 (a = 0, 1, 2, 3); the Hermitian 75 matrix 

75 = 7 5 = 7°7 1 7 2 7 3 = 70717273, (3) 
7! = 75 and 75 = 1; the Hermitian spinor matrix, 

<? ab = l -[i a ,i b ]- (4) 



Totally antisymmetric tensor e pupa = e abc de p e b e p c e d . The space-time metric of four-dimensional 
Euclidean manifold 1Z A is 

9iiv (x) = e p {x)e b {x)5 ab = ~^{e p ,e u }. (5) 

And the Lorentz scalar components of the metric tensor are then simply 

S a b = ^ a e y h , (6) 

where the inverse of the tetrad fields e^ a e v a = 5^ v and e b e^ a = 5 b a . 

Two gauge invariances due to the equivalence principle have to be respected: (1) the diffeo- 
morphism invariance under the general coordinate transformation x — > x'(x); (2) the local gauge 
invariance under the local Lorentz coordinate transformation — > £,'(x), i.e., 

C' a (x) = [A(x)]tt b (x)- (7) 
Under the local Lorentz coordinate transformation (|7|), the finite local gauge transformation is 

V(£)=expi[0 a6 (£K 6 ] G <?0(4) 

V(07aVt(0 = [A- 1 (x)]^ b , (8) 



4 



where 9 ab (^) is the antisymmetric tensor and an arbitrary function of £ = £(x). The Dirac-matrix 
valued fields e^, oj^ and fermion spinor field ip are transformed as follows 

e„(0 = v(0^(Ov f (0; (9) 

-»• ^(0 = v(0«A.(0v f (0 + v(0 Vt(£), (io) 

V>(£)^V>'(6 = v(0V-(6; (ii) 

p; = v(o^v f (0, (12) 

where the derivative = e^(9/3£ a ), the covariant derivative 

V^ = d^-iguj^0, (13) 

and g is the gauge coupling. Corresponding to the finite local gauge transformations (|9lllip. in- 
finitesimal local gauge transformations for fields e^, and ip are 

MO = ab (646, c e£(6; (14) 

Scj^) = 2 75 e abcd ojf9 cd (0-io- ab d^ ab (0; (15) 

WO = iO a \0^(0, (16) 

where 

dab,c = i[(7ab, 7c] = 2(<5fc c 7a - (5ac7fe)> (17) 

and we use the commutator relation 

{a a P,a^} = -2ij 5 e a ^, (18) 



to obtain Eq. 0151) . 

In an SU(2) gauge theory, gauge field -A (£e) can be viewed as a connection J A a (^E)d^% on 
the global flat manifold. On a locally flat manifold, the spin connection uo^dx^" = w a (£)<i£ a , where 
w a (£) = Ufi^ai one can identify that the spin-connection field uj^x) or £j a (£) is the gravity analog 
of gauge field and its local curvature is given by 

R ab = du ab -guj ae Auj b e , (19) 

and the Dirac-matrix valued curvature R^ = RfycTab- Under the gauge transformation (|9ll0p . 

R ab = V(0R ab (0V\0- (20) 
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The diffeomorphism invariance under the general coordinate transformation x —> x'{x) is preserved 
by all derivatives and d-form fields on TZ A made to be coordinate scalars with the help of tetrad 
fields e^ a = d^/dx^ (see Ref. [llj]). The diffeomorphism and local gauge- invariant EC action for 
gravity coupling to fermions is given by the Palatini action Sp and host modification Sh for the 
gravitational field, 

S E c(e,uj) = S P (e,uj) + S H (e,uj) + S F (e,uj,ifj) (21) 
S P {e, u) = j-J d*x det{e)e abcd e a A e b A R cd , (22) 

S H (e,u) = — [d 4 xdet(e)e a Ae b AR ab , (23) 



and fermion action Sf (see Refs. [Hj 



S F (e,u,i/>) = ^Jd 4 xdet(e)[ije fl V^ + h.c.], (24) 

where k = 87rG, the Newton constant G = l/wpj anck , det(e) is the Jacobi of mapping x — >• 
and the integration J d 4 x = J^ 4 d 4 x. The complex Ashtekar connection [fj]] with reality condition 
and the real Barbero connection UM are linked by a canonical transformation of the connection 

n n 

with a finite complex Immirzi parameter 7 ^ [15l |. which is crucial for loop quantum gravity [16]. 

Classical equations of motion can be obtained by the stationarity of the EC action ([2T|) under 
variations (IMTTTl. 



5S E c(e,uj,ip) = : EC Sen + S f° Sip(x) + SeC 6u)n = 0. (25) 
de^ r dip(x) dw M 

From Eqs. (11411160 . we find that Eq. (|25[) can be expressed in terms of independent bases 75, 7^ 

and a a b of Dirac matrices. Therefore, for arbitrary function 9 a b(£,)> Eq. (1250 leads to the following 

three equalities 

^ = 0; ^ = 0; ^ = 0. (26) 
The first and second equations respectively lead to the Dirac equation, 

e"V^(x) = 0, (27) 

and the Einstein equation 

eabcde a A e b A R cd [uo{e)} = /#(x)(e A V)i>(x), (28) 
where the energy-momentum tensor is 

${e A V)ip = -^V v - V^ty. (29) 
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The gauge invariance of the EC action (|2ip under the gauge transformation (|15p leads to the third 
constraint equation SSec/^h = °f Eq. (j26l) . which is the Cartan structure equation, 

de a - guj ab A e 6 - T a = 0, (30) 

where the nonvanishing torsion field, 

T a = K ge b A e c J ab ' c , (31) 

relating to the fermion spin current 

= i^{a ab , 7 C }V> = e abcd ^ d7 5 i>, (32) 

W ab n c } = ie ahcd ^ ld . (33) 

The fermion spin-current (|32|) contributes only to the pseudotrace axial vector of torsion tensor, 
which is one of irreducible parts of torsion tensor [171 ]. The solution to Eq. (130ft is 

< = <(e) + <, < = K 5 e^J afe c , (34) 

where the connection io ab (e) obeys Eq. ()30|) for torsion- free case T a = 0, 

de a -guj ab (e) Ae b = 0. (35) 

Replacing the spin-connection field in the Einstein-Cartan action (|22I24|) by Eq. (|34p . 

5p[e,w] Sp[e,w(e)] + ng 2 J d^de^e)^-/ 5 ^)^^ 5 ^); (36) 
S F [e,w,^] -> 5 F [e,^(e),V',^ + 2^ 2 y'd 4 xdet( e )(^7VV')(^7 5 ^), (37) 

one obtains the well-known Einstein-Cartan theory: the standard tetrad action of torsion-free 
gravity coupling to fermions with four-fermion interactions, 

SEc[e,u(e),ip,'ip] = S P [e, w(e)] + S F [e, w(e), ip, if] 

+ 3ng 2 f d 4 xdet(e)(^7VV')(V'7d7 5 V')- (38) 



Note that the four-fermion interaction actually is the coupling of two fermion spin-currents (J32J). 
Taking into account the host action (|23|) . one obtains 

S E c[e,u(e),ip] = Sp[e,u>(e)] + S H [e,u>(e)] + S F [e,uj(e),ip] + S iF (e,tp) (39) 
S 4F (e,ijj) = K^g 2 f d 4 xdet(e)(^ 7 VVO(VWV), (40) 



where £ = 7 2 /(7 2 + 1) [3]. Using the commutator relations ([T5]) and [er a &j75] = 0, one can show- 
that (V'7d7 5 V') is a pseudovector and (I40D is invariant under the gauge transformation (jilt) . 

As we can see from Eqs. (|24p to (]39p . the bilinear term (I24p of massless fermion fields coupled 
to the spin-connection field (|13p is bound to yield a nonvanishing torsion field T a (|30p . which is 



local and static (see, for example, Refs. 12j, [l9|]). As a result, the spin connection is no longer 



torsion- free and acquires a torsion-related spin-connection Gj^ (I34p . in addition to the torsion- 
free spin connection w^ 6 (e). The torsion-related spin connection oof is related to the fermion 
spin current (|32p . The quadratic term of the spin connection field uj in the curvature ()19p and 
the coupling between the spin connection field u and fermion spin current in Eqs. ()13|24p lead 
to the quadralinear terms of fermion fields in Eqs. (|36p and (|37p . Another way to see this is to 
treat the static torsion-related spin-connection uf ([HI) as a static auxiliary field, which has its 
quadratic term and linear coupling to the spin current of fermion fields. Performing the Gaussian 
integral of the static auxiliary field, we exactly obtain the quadralinear term (|40p . in addition to 
the torsion-free EC action. 

The action (j2Tj) and classical Eqs. (|27ti30p can be separated into left- and right-handed parts 
[20I . with respect to the local SUl(2) and SUr(2) symmetries of the Lorentz group £0(4) = 
SUl(2) ® SUr(2). This can be shown by writing Dirac fermions ip = ipL + ipRi where Weyl 
fermions tpL,R = Pl,riP, Pl,r = (1 T 7s)/2; and Dirac-matrix valued tetrad field = + e R , 

R = Pl^r&^i as well as Dirac-matrix valued spin-connection fields uj^ = uij^+ui R ,ui^ R = P^w^. 



III. THE REGULARIZED EINSTEIN-CARTAN THEORY 



A. Simplicial complex 

The four-dimensional Euclidean manifold 1Z 4 is discretized as an ensemble of A/q space-time 
points (vertexes) "x E and Ai links (edges) "/^(x)" connecting two neighboring vertexes. 
This ensemble forms a simplicial manifold M embedded into the 1Z 4 . The way to construct a 
simplicial manifold depends also on the assumed topology of the manifold, which gives geometric 
constrains on the numbers of subsimplices (Ao, Ai, • • •, see Ref. Qj). In this article, analogously to 
the simplicial manifold adopted by the Regge calculus we consider the simplicial manifold A4 as a 
simplicial complex, whose elementary building block is a 4-simplex (pentachoron). The 4-simplex 
has five vertexes - 0-simplex (a space-time point "2"), five "faces" - 3-simplex (a tetrahedron), and 
each 3-simplex has four faces - 2-simplex [a triangle h(x)], and each 2-simplex has three faces - 
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Ul(x + a u ) 
e^{x + a u ) 




/ e v (x + ) , U v (x + ) 



"/ x + a M + a 



V 



FIG. 1: We sketch a 2-simplex (triangle) h(x) formed by three edges Z M (x) = ae jU (x), l p (x + a p ) — ae p (x + a p ) 
and Zj,(x + a ;i ) = ae v {x + a„) [a = 1] connecting three vertexes x, x + a p and x + a„. Assuming three 
edge spacings a pi a„ and a p (|4ip are so small that the geometry of the interior of each 4-simplex and its 
subsimplex (3- and 2-simplex) is approximately flat, we assign a local Lorentz frame to each 4-simplex. On 
the local Lorentz manifold £ a (x) at a space-time point "x", we sketch a closed parallelogram Cp(x) lying 
in the 2-simplex h{x). Its two edges e p (x) and e£,(x) are two edges of the 2-simplex h(x), and other two 
edges (dashed lines) e^(x + a„) and e„(x + a M ) are parallel transports of ej^x) and ej,(a;) along v and /Li 
directions, respectively [see Eqs. (I46I47P and (|62I63I) ]. Each 2-simplex in the simplicial complex has a closed 
parallelogram lying in it. Group- valued gauge fields U p (x) and U^(x) are respectively associated to edges 
e /J (x) and ej,(x) of the 2-simplex h(x), as indicated. The fields e p = e p (x + a p ) and U p = U p (x + a^) are 
associated to the third edge (x + a M , p) of the 2-simplex h(x). The group fields U u {x + a M ) and Ut(x + a v ) 
indicate the parallel transports of U^(x) and U^x) [see Eqs. ()48I49[) and (I82I83P ] for the zero curvature 
case. Note that the point (x + a p + a v ) is not a vertex of the simplicial complex, points: {x — a p ) 1 {x — a„), 
(x + a p + a p ), (x + a M — a p ), and (x + a v + a p ), which are not shown in the sketch, are not vertexes of 
the simplicial complex as well. Parallel transports e„(x + a p ) and e) p (x + a v ), as well as U„(x + a p ) and 
Ul(x + a v ) are not associated to any edge of the simplicial complex. Throughout this article, the notations 
e and U indicates parallel transports that are not associated to any edge of the simplicial complex. 

1- simplex [an edge or a link "Z„(aj)"]. Different configurations of the simplicial complex correspond 
to variations of relative vertex-positions {x}, edges u {lu(x)} n and "deficit angles" associating to 

2- simplices h(x). These configurations will be described by the configurations of dynamical tetrad 
fields e At (x) and spin-connection fields uJ^ix) assigned to 1-simplexes (edges) of the simplicial com- 
plex in this article. We are not clear now how to relate configurations of fields e^(x) and oj^x) to 
topological constrained configurations of the simplicial complex in dynamical triangulations. 
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1. Edges: 1-simplexes 

The edge (1-simplex) denoted by (x, connecting two neighboring vertexes labeled by x and 
x + a M , can be represented as a four- vector field l^(x), defined at the vertex "x" by its forward 
direction \x pointing from x to x + a M and its length 

a M (x) = |^(x)|^0, (41) 

which is the distance between two vertexes x and x + a^. The fundamental tetrad field e«(x) is 
assigned to each edge (1-simplex) of the simplicial complex to describe the edge location "x," di- 
rection "/u" and length a M (x). We use the tetrad field e M (x), defined at the vertex x, to characterize 
the edge (1-simplex) l^{x) 

l^x) = ae M (x), (42) 

where the Planck length a = (SttG) 1 ^ 2 = k 1 / 2 , and 

Ux)\ = ^{\tT[ efl (x)- efl (x)]\} 1/2 . (43) 

By definition, either l^x) or e^(x) is a Dirac-matrix valued four- vector field, defined at the vertex 
V. 

2. Triangles: 2-simplexes 

We consider an orienting 2-simplex (triangle) (see Fig. [I]). This 2-simplex (triangle) has three 
edges connecting three neighboring vertexes that are labeled by x, x + a M and x + a u . This triangle 
(2-simplex) has two orientations: (i) the anti-clocklike h(x) [ x i — > x + h-^-> a; + a v h^-> x ] and 
(ii) the clocklike h\x) [x i — ^ a; + i — %■ x + i — ^ x]. 

Along the triangle path of the anti-clocklike 2-simplex h{x) [x x + x + a v x], 

three edges and their forward directions are represented by: (1) lu{x) and [i pointing from x to 
x + a M ; (2) l p (x + a M ) and p pointing from x + a M to x + a u ; (3) l u (x + a v ) and ^ pointing from 
x + a u to x. The lengths of three edges are respectively represented by edge spacings o^, a p and 
a u [see Eqs. (|41I43|) ] . We use the tetrad fields 

e M (x), e p (x + a M ), e u (x + a u ), (44) 

defined at x, x+a^ and x + a^, to respectively characterize locations, forward directions and lengths 
of three edges: ([42]) and 

Ipix + an) = ae p (x + a M ), 
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l v (x + a„) = ae u (x + a u ), (45) 
of the anti-clocklike 2-simplex h(x) [see Fig. Q] and Eqs. (I42( I43j) ] . 

B. Parallel transports and curvature 

The fundamental spin-connection fields {oj^(x)} are assigned to 1-simplices (edges) of the sim- 
plicial complex, i.e., each edge (x,/u) we associate with it uj^x). The torsion-free Cartan Eq. (|35l) 
is actually an equation for infinitesimal parallel transports of tetrad fields e"(x). Applying this 
equation to the 2-simplex h(x), as shown in Fig. [IJ we show that e"(x) [e"(x)] undergoes its parallel 
transport to e"(x + a M ) [e"(x + a„)] along the /i (y) direction for an edge spacing a fl (x) [a u (x)], 
following the discretized Cartan equations 

e a v {x + o M ) - e a u (x) - a^guf{x) A e vb {x) = 0, (46) 
^ l {x + a v )-el{x)-a v gojf{x)he tlh {x) = 0. (47) 

The parallel transports e^x + a^) and e®(x+a u ) are neither independent fields, nor assigned to any 
edges of the simplicial complex. They are related to ej,(x)[e^(x)] and uj^{x)[ui u (x)] fields assigned 
to the edges (x, —v) and (x, fj,) of the 2-simplex h{x) by the Cartan Eq. (|46l47p . Because of torsion- 
free, e^(x),el(x) and their parallel transports e^(x + a v ),e v [x + a M ) form a closed parallelogram 



C P (x) (Fig. QJ. Otherwise t 
locally by a flat space-time 



lis would means the curved space-time could not be approximated 
211 ] . Note that the point (x + + a^) at the closed parallelogram 



Cp(x) (Fig. [T]) is not any vertex of the simplicial complex. 

For the zero curvature case R^{x) = 0, the curvature Eq. (|19p can be discretized as, 

u;f(x + a,) - uf(x) - a,g^(x) A J ev {x) = 0, (48) 
ojf(x + a v ) - uf(x) - a v gu?(x) A u%{x) = 0, (49) 

where uj„ b (x + a^) and 0^ b (x + a v ) are respectively parallel transports of ojy b (x) and ui^(x) in the 
fj,- and ^-directions. Analogously to the parallel transports e°(x + a^) and e^(x + a^) given by 
Eqs. (j4"6|) and (jUJ), parallel transports co® b (x + a^) and c^f(x + a v ) are neither independent fields, 
nor assigned to any edge of the simplicial complex. They are related to o^(x) and co u (x) fields 
assigned to the edges (x, fj,) and (x + a u , u) of the 2-simplex h{x) by the parallel transport Eqs. (j4~8|) 
and (|4"9j) . The fields uj^(x), uj u {x) and their parallel transports ^(x + a u ),u) u (x + a M ) also form 
a closed parallelogram, analogously to the one Cp(x) formed by the tetrad fields e^(x),e u (x) and 
their parallel transports e^(x + a u ), e u (x + a M ) (see Fig. [1]). 
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Whereas, for the nonzero curvature case Rf^(x) / 0, the curvature Eq. JT^ can be discretized 

as 

u>f(x + a^-L } f(x)-a^(x)Atvl(x)=a^R*l(x), (50) 
tf(x+a v )-tf(x)-a v gu,?(x) A <(*) = a„i^(s), (51) 

which define fields u® b (x + a M ) and (a; + a„) in terms of fields uj® b (x), u)® b (x) and curvature 
R® Jx). These fields (x+a M ) and (aH-a„) are neither independent fields, nor assigned to any 
edge of the simplicial complex. They are related not only to uj^ b (x) and ui% b (x) fields assigned to 
the edges (x, //) and (x + a u , v) of the 2-simplex h(x), but also to the curvature R ab u (|50l) and i?^ 

(EU). 

These fields u® b (x + a^) and uj^ b (x + a„) are no longer parallel transports Cu® b (x + a^) and 
u>^ b (x + a u ) defined by Eqs. (f48l) and (f49l) . The difference between u^ b (x + a M ) and u% b (x + a M ) [or 
between io ab {x + a v ) and co^ b (x + o„)] is the curvature apR°~ v (x) [a v R®Jx)], 

w^(x+a M )-w?(x+a M ) = afX R%{x), (52) 
<(x+a,) -^(x+a,) = a I/ fl$(s). (53) 

The fields uj^x), uj u (x) and fields w M (x + a u ),uj u (x + a^) do not form a closed parallelogram, due 
to the nonzero curvature R^ b (x) ^ 0. 



C. Group- valued fields 

Instead of a to^x) field, we assign a group-valued field U^x) to each edge (1-simplex) of the 
simplicial complex. On the edge (x, [/,) connecting two vertexes x and x+a^ in the forward direction 
fjL, we place an 50(4) group- valued spin-connection fields, 

U^x) = e i9 ^ (x \ (54) 

whereas the same edge (x + a M , — p) in the backward direction —[i, we associate with it 

U^(x + a^ = Ul(x) = U-\x), (55) 

analogously to the definition of link fields in lattice gauge theories. On the three edges in forward 
directions (x,p), (x + a^, p) and (x + a u ,v) of the anti-clocklike 2-simplex h(x) (/j, ^ v ^ p see 
Fig. [I]), we define 50(4) group- valued spin-connection fields, 

U^x) = e i9aw ^ x \ (56) 
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U p (x + a p ) = e W(x+a M ) ? (57) 

U v {x + aJ) = e i90Wv{ - x+a "\ (58) 

which take values of the fundamental representation of the compact group SO {A). On the three 
edges in backward directions (x, —v), (x + a u , —p) and (x + a p , —fi) of the clocklike 2-simplex fv(x) 
(see Fig. [TJ, we define 50(4) group- valued spin-connection fields, 

U- V {x) = Ul{x + a u ) = e- i9a ^ x+a "\ (59) 
U. p {x + a u ) = U}(x + a p ) = e -^(-+^), ( 60 ) 
C/_ M (x + a M ) = C/t(x) = e -«W. (61) 

These uniquely define group-valued spin-connection fields on the anti-clocklike and clocklike 2- 
simplex. 



1. Unitary operators for parallel transports o/e M (x) fields 

Actually, these group-valued fields (I56H58P and (|59ti6ip can be viewed as unitary operators for 
finite parallel transportations. The parallel transportation (Cartan) Eqs. (146p and (147ft can be 
generalized to (/j, ^ v) 

e u (x + a p ) = Ul(x)e v (x)U^(x), (62) 
e p (x + a v ) = Ul(x)e fl (x)U u (x), (63) 

and using Eq. ([55]) these equations can be equivalently rewritten as 

e„(x) = UL^{x + a p )e v {x + a p )U- p ,(x-\-a p ), (64) 
e p (x) = UL v (x + a v )e p (x + a v )U- v (x + a v ). (65) 

While for (fj, = u), we similarly have the following parallel transportation equations 

e p (x + a p ) = Uj l (x)e p (x)U p (x), 

e p (x) = UL^x + a At )e At (x + a p )U- p (x + a p ), (66) 

indicating that e p (x) is parallel transported to e p (x + a p ) in the fi forward direction, and e p (x + a p ) 
is parallel transported to e p {x) in the —fi backward direction. Similar discussions can be made for 
parallel transports with the unitary operator U p (x + a p ). 
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2. Unitary operators for parallel transports of el(x) fields 

In the simplicial complex, each edge (1-simplex) connecting two vertexes has only one direction. 
One can identify each edge by its starting vertex and direction pointing to its ending vertex. On 
the basis of the tetrad field e At (x) (jl2j) defined at the vertex "x" for the edge (x, fi) starting from the 
vertex "x" in the forward direction {jj) to the vertex "x + a^," below, using the unitary operator 
Ufj,(x) for parallel transports, we will uniquely introduce the "conjugated" field ej^(x) defined at 
the vertex "x" to describe the same edge (x + a M , — fi) but in the backward direction — fj, starting 
from the vertex "x + a^" to the vertex "x." Analogously to Eq. (|42p . this edge starting from the 
vertex "x + a^" in the backward direction (— fi) can be formally represented by 

l-fj,(x + a M ) = ae_ M (x + a^). (67) 

By the parallel transport, we define the field e_ M (x + a^) as 

e_^(x + a^) = Ul(x)ej l (x)U lx (x) = el l (x + a fJ ,), (68) 

in terms of the unitary operator U^x) and conjugated tetrad fields ej^x) defined at the vertex 
"x". From the definition in Eq. (I68p . we rewrite 

ej,(x) = U ll (x)e- fl (x + a fl )Ul(x) = e- fl (x). (69) 

The second equalities in Eq. (|68p and Eq. (I69p are given by the definition of parallel transports by 
unitary operators [see Eq. ([62]) ]. Equation ([68]) means that we can associate the conjugated field 

e+(x) = t/ M (x)4(x + a^Ufe), (70) 

with the same edge (x + a^, — ^) but in backward direction — /i and write 

ll(x) = ael(x). (71) 

As a result, the edge (x,fi) [(x + a^, — /i)] in the forward (backward) direction is uniquely described 
by the field e M (x) [ej^x)] defined at the vertex x. Note that the conjugated field e^(x) is given by 
the parallel transport (jTUl) from x + to x in the direction (— fi). In addition, Eqs. (j68j) and (|69l) 
indicate that conjugated fields mean the inverse of field's direction (/u — > — /u). 

This prescription shows that the edge (x,/i) is completely described by the fields e M (x) and 
e|j(x), latter is a function of fields e M (x) and U^(x), as required by the principle of local gauge 
symmetries and the gauge field U^{x) corresponds a parallel transport between x and x + a^. In 
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consequence, any edge (1-simplex) of the simplicial complex is uniquely identified by its location 
and direction (z,cr), and described by the fields e a (z) and U a (z). 

Using the properties (70)^ = — 7a [see Eq. ([2])] and the definition of tetrad field e /1 (x) = e^ L (x) , y a , 
where the index p is fixed, we have 

et(x) = [e/(x) 7 J t = (7a) t [<(x)]t, 

= -e^x), (72) 

where because of the index p being fixed, the real tetrad-field component e^(x) = d^ a /dx^ can be 
viewed as a one-row matrix (e^ , e^ 1 , e^ 2 , e^ 3 ) and [e^ a (x)]^ a one-column matrix (e^ , e^ 1 , e^ 2 , e^y . 
Analogously to Eq. ([4~3]) . the length of the edge ([71]) in backward direction — p,, 



|lj0«0l = ^{|tr[et(x)-et(x)]|} ] 



2 

= IWI- ( 73 ) 

which is the same as the length of the edge in the forward direction p. 

We turn to the discussion of other two backward-direction edges (x + a u , — v) and (x + a M , — p) 
of the clocklike 2-simplex h){x) (see Fig. [[]). Analogously to Eqs. (|68|69p . we have in the (— u) 
direction, 

e^ v {x) = U v (x)4(x + aM{x) = el(x), (74) 
el(x + a u ) = Ul(x)e- v (x)U v (x) = e- v (x + a u ), 

and in the (— p) direction 

e- p (x + a u ) = U^x + dp)ep(x + ap)U p (x + cip) = e p (x + a u ). (75) 
e\{x + afj) = U p (x + a fJ )e^ p (x + a u )Uj ) (x + a ll ) = e^ p (x + a tJi ), 

As a result, the edge (x + a u ,u) [(x + a^, — u)\ in the forward (backward) direction is uniquely 
described by the field e v {x + a v ) [el(x + a u )] defined at the vertex x + a u 

= Ut(x)el{x)U v {x), (76) 

see Eq. (}T4"|) . Note that the conjugated field e\,{x + a v ) is given by the parallel transport ([76]) from 
x to x + a u in the direction (1/). We can write 

ll(x + a u ) = ael(x + a u ). (77) 
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Similarly, the edge (x+a M , p) [(x+a^, —p)] in the forward (backward) direction is uniquely described 
by the field e p (x + a p ) [et(x + a p )] defined at the vertex x + a p 

e p {x + o M ) = U p (x + a A1 )eJ ) (x + a u )Uj(x + a M ), (78) 

see Eq. ([75]) . Note that the conjugated field e p (x + o M ) is given by the parallel transport ([78]) from 
x + a u to re + a p in the direction (— p). We can write 

l P {x + a p ) = ae p {x + a p ). (79) 

This prescription shows that the edge (x + a u , v) is completely described by the fields e u (x + a u ) 
and U v (x + a v ), and the edge (x + a^p) by the fields e p (x + a p ) and U p (x+a IJi ). The field U u (x + a u ) 
[U p (x + a M )] corresponds a parallel transport between x and x + a p (x + a p and x + a w ). 

Along the triangle path of the clocklike 2-simplex h)[x) [x i — x + a v i — -> x + a p i — ^> x] (see 
Fig. [1]), these three edges and their backward directions are formally represented by (1) Z_ p (x + a p ) 
and — p pointing from x + a /t to x; (2) l~ u {x) and — v pointing from x to x + a u ; (3) Z_ p (x + a u ) 
and — p pointing from x + a v to x + a /t . Based on Eqs. ([55]) . ([75]) . ([75]) . ([TP]) . ([76]) and ([75]) . we use 
the conjugated tetrad fields 

(x + gv), eT(x + a M ), (80) 

which are respectively defined at vertexes x, x+a„, x+a^, to characterize both backward directions 
and lengths of three edges ([7T]) . ([77]) . and ([79]) of the clocklike 2-simplex }v{x). 

In the simplicial complex, each edge (1-simplex), described by tetrad field e fl (x), is uniquely 
identified by its location and direction (x,/x), and each triangle (2-simplex) h(x) has a definite 
orientation, as indicated in Fig.[T] either ant i- clocklike or clocklike. Thus each triangle, for example, 
the one presented in Fig. CD is completely described by the tetrad fields e M (x), e u (x + a u ), e p (x + a p ), 
and unitary operators C/^(x), U u (x + a u ), U p (x + a p ). 

3. Unitary operators and curvature 

In the zero curvature case, the group-valued fields for parallel transports w p (x + a u ) and uj u (x + 
a^), defined by parallel transport Eqs. ([55]) and ([59"]) . are given by 

U p (x + a u ) = e w(^) > u v {x + a u ) = e fcaM*+«iO. ( 8 l) 
Similarly to Eqs. (I62|63p . the parallel transport Eqs. ([48ft and (I49p can be generalized to 

^(x + a M ) = Ul(x)U„(x)Up,(x), (82) 

^(x + a,) = C^(s)?7 /x (x)C^(x). (83) 
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The parallel transport fields U v [x + a^) and U^x + a u ) together with U^{x) and U u (x) form a 
closed parallelogram, see Fig. [TJ This closed parallelogram is not the same as the parallelogram 
Cp{x) formed by e and e fields. 

In the nonzero curvature case, corresponding to the fields uj^(x + a u ) and ui u (x + a fJ ) defined by 
Eqs. (1501) and (|5T1) . the group- valued fields can be similarly given by 

U^X + a v ) = e i9^(x+a v )^ Uu ( x + a j = e i gaUv (x+a lt ) ^ ^ 

whose values obviously depend on the curvature R^ u (x). The same as the fields uj^x + a u ) and 
oj v (x + a^), these group-valued fields U v (x + a M ) and U^x + a u ) are neither independent fields, 
nor assigned to any edge of the simplicial complex. They are related to U^x) and U u (x) fields 
assigned to the edges (x, fi) and (x, u) of the 2-simplex h(x) by 

U v {x + a^) = Ul(x)U u (x)U^x), (85) 
U^x + au) = Ul(x)U^(x)U u (x), (86) 

which are generalized from Eqs. (|50|) and (|5T1) . The fields U u {x + a^) and U^(x + a u ) defined 
in Eqs. (|85l86p encode the information of a nontrivial curvature. They do not form a closed 
parallelogram together with U^x) and U u (x), at the point (x + a M + a u ) (see Fig. [[]). 

In order to see the nontrivial curvature information encoded in the fields U u (x + a^) and U^(x + 
a u ) defined by Eqs. (f8"H) - ([86"j) . based on Eqs. ([85]) and (j86|) . we introduce quantities 

U^ix) = U^U^x) = U^(x)U u (x + o M ), (87) 
^(s) = Up{x)U v (x) = U u (x)U„(x + O, (88) 

and calculate their expressions in the naive continuum limit. In the naive continuum limit: aguo^ <C 
1 (small coupling g or weak uj^ field), indicating that the wavelengths of weak and slow- varying 
fields uj^{x) are much larger than the edge spacing a^, we obtain (see Appendix [A]) 

U^ u {x) = exp ^iga[uj^(x) + u v {x)\ + iga 2 d^ u (x) 

- i( 5 a) 2 [^ t (x),^(x)] + 0(a 3 )}, (89) 

where C(a 3 ) indicates high-order powers of agoj^. It is shown that the quantity U liv {x) [Eq. (I89p ] 
is related to the curvature R^ u (x) in Appendix [Aj For the sake of simplicity in the following 
calculations to show the naive continuum limit, the quantities introduced by ([8"7|) and ([H5|) and 
their expressions in the naive continuum limit (|89|) are useful. 
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D. Triangle constrain and area 

Three tetrad fields e p (x), e p (x + a p ) and e u (x + a u ) [see Eq. (14"4"1) ] are three edges of the anti- 
clocklike 2-simplex h(x), satisfying the triangle constraint 

e p (x + a M ) = e_„(a;) - e^x) = el(x) - e p (x). (90) 

Equivalently, three tetrad fields e 1 (x), e\{x + a v ) and ej,(x + a p ) [see Eqs. ([7D|). ([76]) . and ([78]) or 
(|80p ] of the clocklike 2-simplex h*(x), satisfying the triangle constraint 

e- p {x + a v ) = e p (x) - e_„(a?) = e p (x) - 4(x), (91) 

where e- p (x + a v ) = e p (x + a v ) [see Eq. ([75]) ]. Also, Eq. (j71]l is used for e_„(x) = e\{x) in the 
second equality of Eqs. (|9QI9ip . Two of three edges are independent for a given anti-clocklike 
(clocklike) 2-simplex h(x) [h^(x)]. 

However, in Eqs. (190|9ip . vector fields defined at different vertexes are related without being 
parallel transported to the same vertex, thus these relationships are not proper and does not 
properly transform under local gauge transformations. This is an exactly essential point of local 
gauge symmetries, that gauge fields U for parallel transports are needed to relate variations of 
gauge freedom at different coordinate points. Using the parallel transport by the unitary operator 
Ua(x), we rewrite the triangle constraint (|90p for the anti-clocklike 2-simplex h(x) as 

U^{x)e p {x + a p )Ul(x) = el(x) - e p (x), (92) 

where in the left-handed side, e p (x + a^) is parallel transported from the vertex x + to the 
vertex x to be related to ej)(x) and e p (x) at the same vertex x in the right-handed side. Using 
e p {x) = U p (x)e p (x + a^lfUx), we rewrite Eq. (j92|) as 

e„(x) + e^(x) + e p (x) = 0. (93) 

Using the parallel transport by the unitary operator U u (x), we rewrite the triangle constraint (|9"Tj) 
for the clocklike 2-simplex h'(x) as 

U u {x)e\{x + a v )Ul(x) = e p {x) - et(x). (94) 

where in the left-handed side e p (x + a u ) is parallel transported from the vertex x + a u to the vertex 
x to be related to ej)(x) and e M (x) at the same vertex x in the right-handed side. Equation ([MP is 
identical to Eq. ([92]) or Eq. ([93]) . if we consider e'Jx) = U u (x)e p {x + a u )Ul{x) and &' p (x) = —e p {x). 
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The proper parallel transports by unitary operators can shift the triangle constrain to other vertexes 
for example, x + and x + a v . 

We are now in the position of discussing the area of the 2-simplex h{x). We define the funda- 
mental area operator of the anti-clocklike 2-simplex h{x) (see Fig. [TJ 

S%,(x) = a 2 e p (x)Ae„ u {x) (95) 

at the vertex x. In addition, we can also define the following area operators 

S^ix + dfj) = a 2 e p (x + a M ) A e_ p (x + a M ) (96) 

at the vertex x + a M , and 

S~ p (x + a v ) = a 2 e v (x + a v ) A e_ p (x + a u ) (97) 

at the vertex x + a v . Using Eqs. (f68]) . ([73]) . and ([75]) . we rewrite the area operators ([95]) - ([97]) of 
the anti-clocklike 2-simplex h(x) as 

S^(x) = a 2 e M (x)Aet(x), (98) 
Sp^x + a p ) = a 2 e p (x + a M ) A e+(x + a p ), (99) 
5 , yp( 2; + a ^) = a 2 e u (x + a u ) A e p (x + a„). (100) 

In the following, we show that area operators ([98]) - ([100p defined at three vertexes x, x + a p , and 
x + a u are universal up to parallel transports by unitary operators. Using Eqs. ([68]) and ([92]) . we 
obtain 

Sw( x + a tJ = ° 2 U t( x )i e U x ) ~ e »( x )P»( x ) 
A [/t^e^t^x), 

= a 2 Ul{x)[el(x)Ael(x)]U p {x), 

= Ul(x)S%(x)U,(x). (101) 

Analogously, using Eqs. ([71]) and ([M]) . we obtain 

5^ p (x + a !/ ) = a 2 Ul(x)e v {x)U v {x) 

A ?7+(x)[e jU (x)-et(x)]i7 J/ (x) 

= a 2 Ul{x)e u {x) Ae p (x)U v (x) 

= l4(x)S*(x)U v (x). (102) 
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In Eqs. (jlOip and (|102p . we use eMi) = —e fl (x), e M (x) A e At (x) = e\(x) A ej^x) = ej^x) A e fl (x) = 



and the same for (fi — > i/). This shows that the area operators (l98j) . (j99j) and (llOOj) defined at three 
vertexes of the 2-simplex h(x) are universal up to parallel transports. 

Therefore, Eq. ([95]) or (f98|) defines the area operator of the 2-simplex h(x) 

.2 



S^ix) = — e^{x)el{x) - e|,(x)e At (x) 



a 2 °"a6 



e«(x)et(x)-e«(x)ej(x) 



(103) 



up to parallel transports. As consequence, the area of the 2-simplex h(x) is uniquely determined 
by 

l 



S h (x) = \S%(x)\, Sl(x) 



tr 



(104) 



Its uniqueness [independence of the vertexes x, x + and x + a u of the 2-simplex i.e., 



5 h (x) = \S (x)\ = \S Jx + a^l = \S (x + a u )\, 



(105) 



can be shown by using Eqs. (|101ll02p . 

In the same way as Eqs. (|95p - (|97p . we define the area operators of the clocklike 2-simplex h)(x): 



s vixi x ) = q2 e ~u{x) A e M (x) 



(106) 



) = a 2 e.-yjyX + cifj,) A e p (x + a^) 
-<Sj M + a M ) = S^(x + a M ), 
a 2 e- p (x + a„) A e^x + a v ) 
-Si Ax + a v ) = Sl\(x + a v ), 



whose directions are opposite to the counterparts of anti-clocklike 2-simplex h(x). However, the 
area of the clocklike 2-simplex h^(x) is equal to the area (|104p . 

Based on the definition of 2-simplex h{x) area (jlQ4l) . we can define a volume element around 
the vertex "x" 



dV(x) = dV h (x), dV h {x) = Sl(x) 

h(x) 



(107) 



where dVh(x) indicates the volume element contributed from a 2-simplex h(x), and J2h(x) indicates 
the sum over all 2-simplices h(x) that share the same vertex x. This definition of volume element 
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(|107p indicates that a 2-simplex h(x) contributes the volume element at its three vertexes x, 
x + a p and x + a u . 

Before ending this section, we note that using the parallel transports (f68l) . (i74|) . and (f75l) . one 
can obtain parallel transports of area operators (|95I96I97[) of triangles (2-simplexes), 

S^{x + a p ) = Ul{x)Sl v {x)U^{x), 

S pv (x + a u ) = ul{x)S%,{x)U v (x), (108) 

which are consistent with the definitions of unitary operators Un{x) and U v (x) for parallel trans- 
ports (I62|63p of edges (1-simplexes). The notation "S^" instead of S pu in the left-handed side 
of Eqs. (I108P indicates that the parallel transport is not associated to any triangle of the 

simplicial complex. 

E. Local gauge transformations 

In accordance with Eq. (|10p . the bilocal gauge transformations of three U fields (|56p - (|58p of the 
anti-clocklike 2-simplex h(x) are, 

U p (x) -> V(x)U p {x)V\x + a p ), 
U u (x + a u ) -> V{x + a u )U u {x + a u )V ] {x), 

U p {x + ap) -> V(x + a p )U p (x + a p )V\x + ai/ ), (109) 

and their inverses (I59p -(|6ip of the clocklike 2-simplex h)(x) transform as 

Ul(x) V(x + ap)Ul(x)vHx), 
U}(x + a u ) -> V(x)C/J;(x + a ; ,)V t (x + a ! ,), 

*7j(x + a M ) -> V(x + a^)C/ p t (x + a M )V t (x + a / ,). (110) 

In accordance with Eq. ([9]), the tetrad fields e p (x), e u (x + a u ) and ep(a?+o^) for the anti-clocklike 
2-simplex /i(x) transform under local gauge transformations 

e M (a?) -»• e^(x) = V(x)e M (x)V t (x), 
e^x + a^) -» e^(x + a^) = V(x + a u )e u (x + a^)V'(z + a v ), 
e p (x + a p ) ->• e^(x + a M ) = V(x + a p )e p (x + a At )V t (x + a M ), (111) 
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respectively at the vertexes x, x + a u and x + a p where they are defined. Using above local gauge 
transformations (I109p - (|llip . we obtain the following local gauge transformations of the conjugated 



These local gauge transformations (|112p of the conjugated fields at the vertexes x, x + a u and x + a p 
are in the same manner as that given by Eqs. (jllip . This means that each edge (1-simplex) l p (x) 
of the simplicial complex is uniquely described by tetrad fields e p (x) and et(x), that are defined 
at the vertex x, and covariantly transformed under local gauge transformation. 

It is worthwhile to mention that the transformations (|112p are just conjugated transformations 
fjlllf) . and consistent with the following local gauge transformations: 



which follow the transformation rules of Eq. (jllip . 

It is shown that the tetrad fields (j44|) and their conjugated fields ([80]) given by Eqs. (|70l) . (|76j) . 
and (|78p . as well as the triangle constraints (|92l94p . are gauge covariant, and properly transformed 
under local gauge transformations (|109p - (|112p . The length (|4"3j) or (f73|) of edges (1-simplexes) is 
unique and invariant under local gauge transformations fllQ9|) - (11 12f) . 

Under local gauge transformations (|109p - (|112p . the fundamental area operators (|98p - (|100p of 
the anti-clocklike 2-simplex h(x) are gauge covariant and transform 



which are consistent with Eqs. (jlOip . (|102p . (|109p . and (|110p . and their counterparts [see Eq. (j!06|) ] 
of the clocklike 2-simplex h'(x) transform in the same manner. The parallel transports (|108p of 
area operators transform consistently with Eqs. (I109p . (IllOp and (I114p . However, the area (1104|) 
of the 2-simplex h{x) is unique and invariant under local gauge transformations. 




e^ p (x + a p ) ->■ e'_ M (x + a p ) = V(x + a M )e_ M (x + a /x )V t (x + a p ), 

e_„(:c) -> e'_„(a;) = V(x)e_ i/ (x)V t (x), 
e_ p (x + ay) ->• e'_ p (x + a v ) = V(x + a v )e- p (x + a^V^x + a„), 



(113) 



SU X ) "> S^(x) = V(x)S^(x)vHx), 
S^ p (x + a u ) ->■ S^ p (x + a„) = V(x + a v )S^ p (x + a v )\>\x + a^), 
^(x + a^) -> SjJ,(ar + a M ) = V(x + a p )S^(x + (x + a p ), 



(114) 
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It is worthwhile to mention that under local gauge transformation (|109p - (|llip . parallel transport 
fields (I62p and (1630 transform locally 



in accordance with local gauge transformations (11 11 1) for tetrad fields. Therefore, the closed parallel- 
ogram Cp{x) (see Fig. [1]), formed by e^(x),e u (x) and their parallel transports e^(x+a v ), e v (x+a fl ), 
is invariant under local gauge transformation. This is consistent with the torsion-free condition for 
the existence of local Lorentz frames at each points of a curved space-time. 

The prescription of using tetrad fields e a (z) and gauge fields U a (z) for parallel transports to 
describe edges (1-simplexes) and triangles (2-simplexes) of the simplicial complex fully respects 
the principle of local gauge symmetries. Therefore, this prescription is independent of a particular 
vertex z, oriented edge l a {z) and triangle h(z), because of the gauge invariance. The formulation 
of defining tetrad fields e a (z) at one of edge endpoints "z" and direction "a," and each triangle 
has a definite orientation is gauge invariant. 

However, the gauge transformation properties of fields U u (x + a M ) and U^x + a u ) defined by 
Eqs. ([85]) and ([86]) . as well as U^ u (x) and U Ufl (x) introduced by Eqs. (|87|) and ([88]) . are very 
complicate under the bilocal gauge transformations (I109p and (jllQj) . This implies that we could 
not use these fields to construct a gauge-invariant object. We need to study the object of three U 
fields, Ufj_(x), U p (x + a M ) and U u (x + a u ) along a closed triangle path of each 2-simplex h(x) (see 
Fig. [I]), which will be discussed in the next section. 



To illustrate how to construct a gauge-invariantly regularized EC theory describing dynamical 
configurations of the simplicial complex, we consider anti-clocklike 2-simplex (triangle) h(x) and 
clocklike 2-simplex (triangle) h\x) (see Fig. [1] and Fig. [2]). 

For simplifying notations, we henceforth do not explicitly write negative signs — — u, —p to 
indicate the backward directions of edges. In terms of the tetrad fields e^(x) and e u (x) of the 
2-simplex h(x) (see Fig. [[]), we introduce the following vertex fields v^ u (x): 



e^(x + a v ) -)• e^(x + a u ) 
e v (x + a^) e' v (x + a^) 



V(x + a„)e M (x + a iy )V t (x + a u ), 
V(x + a^)e u {x + a^)V^(x + a M ), 



(115) 



F. Regularized EC action 



e^{x) = o- ab e a (x)Ae b (x) 



(116) 
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= if ah [e°(x)el(x) - e a u (x)e b p (x) 

= ^[e^{x)e v {x) - e v {x)e^{x)}, (117) 

which have properties: v pu {x) = —v up (x), ti[vn V {x)\ = and vj lL ,(x) = v vp {x) (see Appendix iBj) . 
Under the local gauge transformation (|9lllip . the vertex fields (|1 16[) and (|117p transform locally 
at a vertex x, 

-> V{x)v pu {x)V ] (x), (118) 

which is transformed in the same manner as area operators (|114p . In addition to the vertex field 
e^u(x) (|117p at the vertex (x), we can define in the same way the vertex fields e PfJj (x + a p ) at the 
vertex (x + a^), and e up {x + a u ) at the vertex (x + a u ) of the anti-clocklike 2-simplex h(x) (see 
Fig. [[]). Actually, the vertex fields e pu (x) (|117p . e pp (x + a M ) and e^ p (x + a u ) are related to the 
fundamental area operators S*,(x) ® , S^(x + a,,) {99} and S* p (x + a„) {TOD]), e.g., 

S%{x) = ia 2 e^(x). (119) 

As discussions for three area operators in Eqs. {95]) - (jl03p . only one of three vertex-fields e pu (x), 
£pn(x + a p ) and e vp (x + a u ) is independent for the anti-clocklike 2-simplex h{x). As for an clocklike 
2-simplex h'(x), vertex fields can be obtained by using the relations ej ll/ (x) = e up {x) and e pu {x) = 

^Ufl {x) ■ 

Using the tetrad fields e p (x) and vertex fields v p , u (x) to construct coordinate and Lorentz scalars 
to preserve the diffeomorphism and local gauge invariance, we define a smallest holonomy field along 
the closed triangle path of the 2-simplex h(x) (see Fig. [I]): 

X h (v,U) = tr[v up (x)U Pi (x)v pp (x + a p )U p (x + a p )v pu (x + a u )U v (x + a v )} 1 (120) 

whose orientation is anti-clocklike, as shown the left graphic in Fig. [2j Considering the clocklike 
orientation, as shown the right graphic in Fig. [2] we have 

x ciockiike^ 5 ^ = tr [ v ^(x)U v {x)v vp {x + a„)U p (x + a u )v pp (x + a p )U p (x + a p )] 

= X h (v,U)\^ v . (121) 

On the other hand, 

= tr [U u (x)v vp (x + a u )U p (x + a u )v pp (x + a p )U p (x + a p )v pu (x)] 
= tr [v pu {x)U v (x)v vp (x + a u )U p (x + a v )v pp (x + a p )U p (x + a p )\ 
= X£ locklik >, [/) (122) 
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where in the second line of equation, we use the properties U$,(x + a v ) = U v (x), Uj,(x + a^) = 
U p {x + a v ), Ul(x) = U^{x + a^) and vj iu (x) = v ull (x). Therefore we have 

X h (v, U) + h.c. = X h (v, U) + X$ ocklike (v, U). (123) 

Equations (1121^ - (1123f) are invariant under gauge transformations (|109p . (|1 10[) . and (|1 18[) . 

Using Eqs. (|120p - (|123p . we are ready to construct the diffeomorphism and local gauge-invariant 
regularized EC action. First we consider the case v fll/ (x) = e liu (x)^: 

A P (e,U) = A E {X h (v,U)+h.c.}, (124) 
9 heM 

where YlheM ^ s ^ ne sum over an 2-simplices h of the simplicial complex. In the naive continuum 
limit: agoj^ <C 1, Eq. (|124p becomes (see Appendix |B|) 

A P (e, Sl(x)e cdab e c A e d A R ab + 0(a 4 ), (125) 

a haM 

where the 2-simplex h(x) contributed volume element S^(x) is given in Eq. (|104p or Eq. ()B17p . 
Based the volume element dV{x) (11071) around the vertex "x" 

E Sl(x) = WdV(x) (126) 

where J2x stands for a sum overall vertexes (0-simplices) of the simplicial complex, and the factor 
1/3 is due to each 2-simplex contributing its area to its three vertexes. The interior of the 4-simplex 
is approximately flat, leading to 

J2 dV ( x ) => I d ^i x ) = I d A xdet[e(x)}. (127) 

x J ^ 

As a result, Eq. (|125p approaches to Sp(e,ui) (|22p with an effective Newton constant 

G eS = \gG, (128) 

and K e s = 87rG e g. The second we consider the 

A H {e,U^) = -i- E [X h (v,U)+h.c], (129) 
9 1 heM 

where the real parameter 7 = [see Eq. (I23p ], Analogously, in the naive continuum limit: 
agu^ <C 1, Eq. (|129p approaches to Sh(z, oj) ([23]) [see Appendix |B] . 

A H (e,UJ = r— r [ d 4 xdet{e(x)]e a Ae b AR ab + (D(a 4 ), (130) 
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X h (v,U) 



Xl(v,U) 



FIG. 2: The smallest holonomy field along a closed triangle path of the 2-simplex h(x): the anti-clocklike 
orientation Xh(v, U) [left]; the clocklike orientation Xj l (v 1 U) [right]. 



with the effective Newton constant n e s = 87rG c ff ([128p . The diffeomorphism and local gauge- 
invariant regularized EC action is then given by 



Aec = A P + A H - 



(131) 



In addition, we can generalize the link field U^(x) to be all irreducible representations j of 
the gauge group SO (4). The regularized EC action (|13ip should be a sum over all irreducible 
representations j, 



Ae G = Y.J. [ A W^ U h) + A j H {e^ Up 



(132) 



where dj is the dimensions of the irreducible representations j and dj = 4 for the fundamental 
representation, which is the dimension of the Dirac spinor space. 



G. Invariant holonomy fields along a large loop 

We consider the following diffeomorphism and local gauge-invariant holonomy fields along a loop 
C on the Euclidean manifold 1Z 4 

X c {v,u)) = Vctiexp i^ig j> v ^ v {x)^ {x)dx v \ , (133) 

where Vc is the path-ordering and "tr" denotes the trace over spinor space. We attempt to 
regularize these holonomy fields f)133|) on the simplicial complex A4. Suppose that an orientating 
loop C passes space-time points (vertexes) xi,X2,xs,- ■ -,xn = x\ and edges connecting between 
neighboring points in the simplicial complex M (see the diagram in the left-hand side of graphic 
equation, Fig. S]). At each point Xi two tetrad fields e M (xj) and e /i /(xj) (fj, ^ fx') respectively 
orientating path incoming to (i — 1 — > i) and outgoing from (i — y i + 1) the point Xi, we have the 
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vertex field v^'(xi) defined by Eqs. (|116p and (|117p . Link fields U^{xi) are defined on edges lying 
in the loop C. Recalling the relationship = UUxi) [see Eqs. (|59H61|) ]. we can write the 

regularization of the holonomy fields (j!33|) as 

X c (v,U) = Pc*tr v^ / (x 1 )U fl >(xi)v fl/u (x2)U u (x 2 ) 
■■■ Vp pl (x i )U f /(x i )v p > a (x i+1 ) 

■ ■ ■ Vx^XN-ijU^XN-i) , (134) 

which preserve diffeomorphism and local gauge invariances. The holonomy fields Xc(e,U) are 
functionals of fields (v, U) and loop C. Consistently with the holonomy fields Xc(e, U) [Eq. (|134p j. 
the holonomy field Xh(e, U) [Eq. (I120p ] is the one with the smallest loop, i.e., the closed path of 
the 2-simplex (triangle) h(x), see Fig. [TJ 



H. Euclidean partition function 

The partition function Zec an d effective action A e ^ c are given by 

Z EC = exp -Af c = J VeVU exp -A EC , (135) 
with the diffeomorphism and local gauge-invariant measure 

fveVU= Y[ f de^dU^x) 5(A) (136) 

where Y\i^{ x )£M indicates the product of overall edges (1-simplices) of the four-dimensional sim- 
plicial complex A4. As already mentioned, the configuration {lu(x) € M} is formulated such that 
each edge l^ix) = ae^x) is defined by giving its coordinate (vertex) x € M. in one of the endpoint 
coordinates x and x + a^, and giving its forward direction \i pointing from x to x + a„. This end- 
point coordinate x and forward direction \i have to be uniquely chosen for each edge lu{%) G -M. 
Beside, on such defined edge l^(x), we place an independent gauge field U^(x) corresponding a 
parallel transport between x and x + a^. The gauge invariant properties, discussed in Sec. IIIIEl 
guarantee that the change of a formulation does not lead to the change in the measure of the 
configuration {lu(x) € M}. In addition, the triangle constraint (|92p and (I93p must be imposed in 
the measure (|136p . symbolically indicated as 5(A), a 5 functional of Eq. (|92|) or Eq. (f93l) . 
In the single edge measure [see Eq. (|136j) ] 

/ de^x)dU^x), (137) 
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dUfj,(x) is the invariant Haar measure of the compact gauge group SO (A) or SUl(2) <g> SUr(2), 
and de^(x) is the measure of the Dirac-matrix valued field e^(x) = J2a e< fj,( x )lai determined by 
the functional measure de^(x) of the bosonic field e^(x). The single edge measure has to be the 
measure over fields only e M (x) and U^x) of the edge in the forward direction fi, because ej^x) 
and UUx) of the edge in the backward direction — \i are related to the fields e At (x) and U^{x) by 
Eqs. (l55j) . (|68|) . (j70j) . and (j72j) so that the single edge measure (I137j) is actually over all degrees of 
fields assigned on the edge. 

It should be mentioned that the measure (|136p is just a lattice form of the standard DeWitt 
functional measure 22| over the continuum degrees, with the integral of the spin-connection field 
u^{x) replaced by the Haar integral over the C/„(x)'s, analytical integration or numerical simulations 
runs overall configuration space of continuum degrees and no gauge fixing is needed. In addition, 
it should be noted that the measure (I136P does not contain parallel transport fields e and U, for 
examples e u (x + a^) and e M (x + a u ) (see Fig. [[]) given by the Cartan Eqs. (|4"6j) and since 
parallel transport fields are not associated to any edges of the four-dimensional simplicial complex. 
This means that the torsion- free Cartan equation has been taken into account. 

In this path-integral quantization formalism, the partition function (|135[) presents all dynamical 
configurations of the simplicial complex, described by the configurations of dynamical fields e^(x) 
and U/j,(x) in the weight of exp — Aec- The effective action A c ^ c (1135(1 contains all one-particle 
irreducible (1PI) functions(operators), i.e., all truncated n-point Green- functions. The vacuum 
expectational values (v.e.v.) of diffeomorphism and local gauge-invariant quantities, for instance 
holonomy fields (I134p . are given by 

(X c (v,U)) = -^l VeVu[x c (v,U)] exp-A E c- (138) 

In the action (1124|129p . Xh(v, U) [Eq. (I120p ] contains the quadratic term of e At (x)-field associated 
to each edge of 2-simplex h(x), the partition function Zec f|135j) and v.e.v. (|138p are not divergent 
for large fluctuating fields, provided the action Aec is positive definite, see discussions below. 
On the other hand, all edge lengths do not vanish [|e M (x)[ ^ 0, see Eqs. (|41|42|) ], .and all simplicial 
triangle inequalities and their higher dimensional analogs should be imposed [3, Q|. Integrating 
spin-connection fields over the Haar measure of compact gauge groups is similar to that in the 
Wilson-lattice QCD, the difference is that the Xh(v, U) (I120p contains three U fields in a 2-simplex 
h, while the Wilson action contains four U fields in a plaquette. Equation (|138p can be calculated 
by numerical Monte-Carlo simulations. We are trying do some numerical Monte-Carlo simulations, 
it will take time so that the results will be published in a separate paper. 
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Before ending this section, we make some discussions on the convergences of the partition 
function (I135|) and v.e.v. (|138|) , Suppose that we first integrate Eqs. (I135p and (|138f) over the 
compact Haar measure of the SO (4) gauge group, roughly speaking, the result gives, in addition 
to a polynomial of tetrad fields e, a combination of both decreasing exponents exp [— ^4( + )(e)] and 
increasing exponents exp [— A^~\e)] as functions of increasing tetrad fields e. From the regularized 
action (I120p . one can find that A^\e) depend on 2-simplex area operators Sh (jl04[) and are the 
sum over all 2-simplexes. A^(e) are either some extremal values of the action Aec f 1 1 3 1 [) with 
respect to group-valued U fields, or those values taken at the boundary points of the compact 
SO(4) gauge group. Clearly, for the case of decreasing exponents exp [— „4.( + )(e)], integrations 
Eqs. (|135p and (|138[) over tetrad fields e are convergent. This is certainly the case for perturbative 
weak U fields, i.e., U ~ 1. While for the case of increasing exponents exp [— A^~\ej\, integrations 
Eqs. (I135P and (I138P over tetrad fields e are divergent. 

To avoid these possible divergences, it is necessary to add into the regularized action Aec (|13ip 
an additional term of another dimensionality: either a curvature squared R 2 term: X^(v, U) + h.c. 
with a new coupling parameter; or a bare cosmological term: A\(e). We consider here an additional 
bare cosmological term Aa to the regularized action Aec (|131 j) : Aec - ► Aec + A\, 



+ h.c. 



_A_ 
4 -"(4!) 

A^det[e£(x)] + h.c. (139) 



Aa ^ = T^^ e ^ Pa T, tr [^(x)e u (x)e p (x)e CT (x) 



where the cosmological parameter A = Aa 2 and A is the bare cosmological constant. The bare 
cosmological term A\(e) is a four-dimensional volume term (sum over all vertexes x), which is 
independent of configurations of group- valued U fields. The exponent exp [— -4a(c)] decreases with 
strong tetrad fields e, large volume configurations. Bare parameters g, 7 and A play an important 
role for convergences of the partition function (|135p and vacuum expectational values (|138p . It 
needs further studies to find the region of bare parameters g, 7 and A for the convergences, and the 
scaling invariant region (g c , 7 C; A c ) for the physically sensible continuum limit, see the discussions 
m the last Sec. ELQ 
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I. Local gauge symmetry 

Analogously to Eq. (f25|) . the local gauge invariance of the partition function (1135|) . i.e., 5Zec = 
under the gauge transformation (|109j) and (|118h . leads to (no summation over index /jl) 

(^e M + ^^ + h.c.)=0. (140) 

Based on Se^ and Sco^ (|14p and (|15p for an arbitrary function 9 a ^(x) and the independent bases of 
Dirac matrices 75, 7 M and a a b, we obtain the "averaged" Cartan Eq. (f35l) for the torsion-free case, 



where we use 

^ = togfa*^ ~ Uf-^l (142) 

for the group- valued field U^(x) = expligaui^x)] ([56]) . The averaged torsion- free Cartan Eq. (|14ip 
actually shows the impossibility of spontaneous breaking of the local gauge symmetry. This should 
not be surprised, since the torsion-free (|30p is a necessary condition to have a local Lorentz frame, 
therefore a local gauge-invariance, as required by the equivalence principle. 

IV. INCLUDING FERMION FIELDS 
A. Bilinear and quadralinear-fermion actions 

Introducing dimensionless fermion field ip'(x) = a 3 ^ 2 ip(x) (drop "prime" henceforth) and using 
the relations 7°(7 a ) t 7° = la, 7 < Va6) t 7° = °ab and 

7^7° = e,; 7 °^7° = (143) 

we consider the following regularized kinetic action of fermion fields, 

A F (e, U,iP) = -J2 [^{x)e»{x)U^x)^{x + a„) 

X, [I 

-i){x + a tl )Ul{x)e tl {x)i){x)\, (144) 



where fermion fields ip{x) and ip(x + a^) are defined at two neighboring points (vertexes) of the 
edge (x,x + a M ), (see Fig. [Q), where fields U^(x) and e^(x) are added to preserve local gauge and 
diffeomorphism invariances, and J2 X is the sum over all edges (1-simplexes) of the simplicial 
complex. 
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Using Eq. (|142p and performing a variation of the regularized fermion action (|144p with respect 
to the spin-connection field uj^x), i.e., <L4p(e, U,ip)/5oj^, we obtain the nonvanishing torsion field 
T a = Kgeb A e c J ah ' c , where the regularized fermion spin current is 

^^e^xh^^x + a,), ^ fixed, (145) 

[see Eq. (|32j)]. Instead of solving regularized Cartan equation and finding an effective the- 
ory, as what is done in the continuum case (125p - (132p . we assume that the U^(x) in Eqs. (I144|) 
and ()145p is the group- valued spin-connection field uj^{e) for the torsion- free case ([35]) . i.e., 
U^(x) = exp[iaguj^(e)]. Thus, the regularization of the effective EC theory ([39]) and (j4"0j) is given 
by Eqs. (1131 1) and (I144p and the regularized four- fermion interaction 

>M^) = SC^E^W^V^W^^ + ^^^^ + a^f/t^T 5 ^^)], (146) 

X, fl 

where (, = 7 2 /(7 2 + 1) = 7 2 /(7 2 + 1) [ see Eq. (|40p] . In the naive continuum limit agw^ <C 
1, the regularized fermion action Af(g,U,iP) (|144p approaches to the continuum fermion action 
5"i?(e, w^, (f24"p . and Eqs. (j!45|) and (|146p . respectively approach to their continuum counterparts 
J ab ' c (|32p and S^p ()40p . The diffeomorphism and local gauge- invariant regularized EC action is 
then given by 

Aec = Ap + Ah + A f + Aaf- (147) 
The partition function Zec and effective action A c ^ c are 

Z EC = exp -Af c = J VeVUVi) exp -A EC , (148) 
with the diffeomorphism and local gauge-invariant measure 

heVUViP = Y[ f de^(x)dU^(x)6(A) ■ ]J f dip (x) dtp (x), (149) 

where dip (x) dip (x) is the measure of Grassmann anticommuting fields. Analogously to Eq. (j!32|) . 
Eqs. (|147p - ([149p can be straightforwardly generalized to include all irreducible representations j of 
the gauge group SO(4) that couple to corresponding spinor states of fermion fields. 

B. Holonomy fields with fermions 

We consider the following diffeomorphism and local gauge-invariant quantities 

X c (e,u,ip) = ip(xi)V exp ^ig J^v flu (x)uj fM (x)dx u ^ ip(x N ), (150) 
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where C stands for an orientating (V) path connecting two vertexes x\ and xn (%i 7^ aJjv) on the 
simplicial complex A4. In Eq. (]150p . Xc(e,u,ip) represents the evolution of the spin of fermion 
fields from the vertex xn to the vertex x\ under the gravitational field influence. Analogously to 
discussions in Sec. pil Gp for the holonomy fields (|133p . we regularize these quantities f| 1 50 [) on the 
simplicial complex as follows: 



which preserves diffeomorphism and local gauge invariances. The graphic representation of 
Xc(e, U, ip) can be found in Fig. [5] (see the diagram in the left-hand side of graphic equation). 



Analogously to the discussions in the continuum EC theory (see the end of Sec. [TTJ) , the regu- 
larized EC action (I147P can be separated into left- and right-handed parts. Fermion fields i/j are 
decomposed into their left- and right-handed Weyl fields: V = V'L + ^R an d ^l,r = Pl,r^-, where 
the chiral projector Pl,r = (1 T 7s)/2 and the commutators [a ab , Pl,r] = and [j a, y b , Pl,r] = 0. 
The 4x4 Dirac spinor space is split into two independent left- and right-handed 2x2 Weyl spinor 
spaces. In the chiral representation of matrices 7° and a ab 



= ■^{xi)V\U^(xi)v tJ ,i l/ (x2)Uy{x2) 

•■• V pp/ (Xi)U p >(xi)Vp< a (xi +1 ) 

■ ■ ■ v Xll {x N )Ul(x N ) t/j{x N ), 



(151) 



C. Chiral gauge symmetries 




(152) 



(153) 



where = k o k and a 1 (i = 1, 2, 3) are the Pauli matrices, we define 7^ R = Ph,Rl a '- 




(154) 



and a' 



Pl,ro" 



ah . 



L,R — 




(155) 
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Using Eq. (|154p . we separate tetrad fields e M into their left- and right-handed fields: e M = + e^, 

R = Pl,r^ ■ Using Eq. (|155|) , we separate spin-connection fields and vertex- fields 
into their left- and right-handed fields: = oj^ + = Pl,r^] and v^ u = v^ v + v R u , 

vfy = Pl.rv^v This splits the Lie algebra of the group SO (4) into two independent Lie algebra 
of sub groups SUl(2)®SUr(2). Therefore, the four-dimensional rotational group SO '(4) is split into 
two commuting and independent groups SU L (2)®SU R (2). The link fields U^x) = U*(x)@Ufa), 
where Ufa) G SU R {2) and Ufa) G SU L {2) respectively. 

The regularized EC theory ()14T[) - (jl49p possesses exact chiral gauge symmetries, as consequences, 
the holonomy fields (|120|) . (|134p . and (|15ip can be split into the left- and right-handed parts: 

X h (e,U) = Xk(e L ,U L ) + X*(e R ,U R )- (156) 
X c (e,U) = X^(e L ,U L ) + X R (e R ,U R ); (157) 
X c (e,U,i/>) = X^{e L ,U L ^ L )+X R {e R ,U R ^ R ), (158) 

where notations in the right-handed side of equations, for instance, X^ (e L , U L , iPl) indicates the 
same function Xc(e,U,i/)) (|15ip with replacements e —> e L , U — > U L and ip — > ipL. The fermion 
action (j!44p and four-fermion interaction (|146p are also separated into the left- and right-handed 
parts: 

A F (e,U,il>) = A F {e L ,U L ,i> L ) + A%(e R ,U R ,ijj R ); (159) 
A AF (U^) = A% F (U L ,ip L ) + A R F (U R ,ip R ). (160) 

The chiral gauge symmetries of the regularized EC theory (I147p - (ll49p are crucial for formulating the 
parity-violating (chiral) gauge symmetries SUl{2) <8> Uy(l), e.g., the standard model for particle 
physics, onto the simplicial complex described by the dynamical tetrad fields c^{x) and group- 
valued spin-connection fields U fl (x). We only discuss the case of Weyl fermions (massless Dirac 
fermions), and the discussions on the case of Majorana fermions are the same, thus not presented 
in this article. 



V. DYNAMICAL EQUATIONS FOR HOLONOMY FIELDS 

Under a local gauge transformation (|9]l- (fTT]l . equivalently ([9]), (1109|) . and (fTTP . the local gauge 

invariance of holonomy fields (X) [Eq. (|138|) ] . i.e., 5(X) = 0, leads to the dynamical equations for 

the holonomy fields X h (fl~20D . X c (gM)) and X c (gSTj) , 

/ sx x v 6Aec x v . ,SX 5A E c x/ , 
Se^ p Se^ p dip dtp 
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B 



FIG. 3: We sketch a graphic representation of the dynamical Eq. (|165p for the smallest holonomy field 
Xh(v,U) (|120j) . The diagram in the left-hand side of the graphic equation indicates the first term in 
Eq. (|165l) . The first and second diagrams in the right-hand side of the graphic equation, respectively 
indicate the third and second terms in Eq. (|165|) . Note that A and A' are the same vertex, so are B and 
B' . We indicate the edge l^, where the local gauge transformation is made. In the right-hand side of the 
graphic equation, the summation over all 2-simplices h(l) associated to this edge l^ is made. 

+iag{XSu lt ) - {X^^-Su^) + h.c. = 0, (161) 
where the index /i is fixed, and for the variation SX/doo^ we use Eq. (|142p and the relationship 

ab 0L > ab OUfi 

Analogously to the analysis in Sec. pilip . we obtain the dynamical equations for the holonomy 
fields X = X h ,X c and X c 

- X 5 -^5e^) + h.c. = 0, (163) 
<^<5V - X S -j^-W) + h.c = 0, (164) 

and 

(X) + (x(ul S ^-))-(x(u^))=0. (165) 

Equation (|165p has the same form as the Dyson- Schwinger equation for the Wilson loops in lattice 
gauge theories. In Figs. 1314151 we show the graphic representations of the dynamical Eqs. (|165p 
for the holonomy fields X h (fT20l) and X c (USD, as well as X c (fTBTD . 

VI. MEAN-FIELD APPROXIMATION 
A. Mean-field approach 



In this section, we try to approximately calculate the partition function (|135p . the vacuum 
expectational values of the 2-simplex area (|104p and the volume element (|107p by using the approach 
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FIG. 4: We sketch a graphic representation of the dynamical Eq. (|165|) for the general holonomy field Xc 
(|134D . The diagram in the left-hand side of the graphic equation indicates the first term in Eq. (|165|) . The 
first and second diagrams in the right-hand side of the graphic equation, respectively, indicate the third and 
second terms in Eq. (|165[) . We indicate the edge 1^, where the local gauge transformation is made. In the 
right-hand side of graphic equation, the summation over all 2-simplices h(l) associated to this edge l p is 
made. 



ip(xi) ^(^i) js< ii(xi) 




FIG. 5: We sketch a graphic representation of the dynamical Eq. (|165|) for the field Xc (|151|) . The diagram 
in the left-hand side of the graphic equation indicates the first term in Eq. (|165[) . The first and second 
diagrams in the right-hand side of the graphic equation respectively indicate the third and second terms in 
Eq. (|165[) . Note that A and A' are the same vertex, so are B and B' . We indicate the edge L, where the 
local gauge transformation is made. We also indicate the fermion field ip^x^) at staring point Xn and the 
fermion field ip(x\) at ending point x\ of the path C. In the right-hand side of the graphic equation, the 
summation over all 2-simplices h(l) associated to this edge l p is made. 

of the mean-field approximation. In the regularized action X^iy , U) (|120p associating to the 2- 
simplex h(x) (Fig. [Q), we replace the vertex fields v pp (x + a p ) and v pu {x + a u ) by assuming a 
nonvanishing mean-field value > 0, 

(166) 

where a, (3 are Dirac spinor indexes. The definition of mean-field value (|166p does not depend on 
whether v pp and v pv contain the matrix 75 or not, due to 75 = 1 and [75, o~ a \\ = 0. The mean-field 
value Mh is independent of any specific vertex, edge and 2-simplex of the simplicial complex. Based 
on the definitions of the 2-simplex area (|104p and the volume element (|10T|) , the mean- field values 



{Ml)5 



ad 



\V pp V pu/ 
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for the 2-simplex area and the volume element are given by 

(Sh(x)) 



a 2 M h , 



{dV(x)) = a A N h M£, 



(167) 



where is the mean value of the number of 2-simplices h{x) that share the same vertex. Note that 
in this preliminary calculations in the mean-field approximation, we do not take into account the 
cosmological term (|139p . since the path integrals are convergent (see below) for positive mean- field 
value M h > 0. 

Based on the mean-field value (|166|) . the smallest holonomy field Xh(v, U) (|12U|) is approximated 
by its mean-field counterpart 



X h (v,U) = tilv^xp^Upix + a^U^x + a^Ml 
xl(v,U) = tv[v^(x)U u (x)U p {x + a^U^x + a p )} 



(168) 
(169) 



where using Eqs. (11211) and (1122ft for fi ^ u we obtain xt{v,U). Note that two of three vertex 
fields v(x) in the Xh(v, U) (|120p . i.e., v pp (x + a M ) and v pv {x + a u ) are replaced by their mean-field 
values Mh, and the 2-simplex h(x) shown in Fig. Q] can also be identified by three different indexes 
H 7^ v 7^ p (no summation over these indexes). Equations (1168P and (1169P depend on U p , and the 
fields (e M , U p ) and (e v , U v ) associated to two edges (x, p) and (x, v) of the 2-simplex (triangle) h(x) 
(see Fig. [1]). Using Eqs. (|168|169p . we define the local mean-field action Ah for the 2-simplex h(x) 



A h = ^[x h (v,U)+xl(v,U) 



+ 



8<7 2 7 



X h {v,U) + x\{v,U) 



tr 



dp{x)T^ (x)e p {x) - e^x)^ Jx)e v {x) 



(170) 



where 



i 

w 

1 

8^ 



75 

i 
2 



Mills 



1 



U v {x)U p {x + a v )UUx) 



+ h.c. 



(171) 



The detailed derivation is given in Appendix [Dj In this mean-field approximation, all 2-simplices 
{h(x)} in the simplicial complex A4 have the same local action ()170p . namely, the single 2-simplex 
mean-field action Ah (|17Up and operator F^ (|17ip are independent of the vertex "x". With the 
local mean-field action (|170p , we define the local mean-field partition function 



VUVeexp — Ah, 



(172) 
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where the local mean-field measure is defined by 

/ VUVe = ( dU^dU u dU p de^de u , (173) 

Jh Jh 

for each 2-simplex h. Thus, the regularized EC action Aec (|131|) is approximated by its mean-field 
counterpart, 

Aec = £ A, (174) 

heM 

which is the sum of the mean-field actions Ah over all 2-simplices h. With the mean-field approx- 
imated action (|174p . we define the mean- field approximated partition function 

Zec= II / VUVeexp-A E C= II 2 h, (175) 
heM h heM 

which is the mean- field counterpart of the partition function (|135[) . 

Using the mean-field EC action Aec (|170p and partition function Zec (|175p . we have the 

following identity 

Zec = Z E c{e- {AEC - AEc) }o, (176) 
where (• • -) is the vacuum expectational value with respect to the mean-field partition function 



Zec (|175p . Using the convexity inequality 



23] 



/ e -(AEC—AEc)\ > e -{AEC-AEc)° ^ (177) 
one can derive the following inequality 

- In Z EC < - In Z EC + (A EC - Aec)o, (178) 

where — In Zec an d — In Zec are proportional to the free energies. We define the right-handed 
side of the inequality (|178p as an approximate free energy (or approximate effective action) 

J^(M h ,g, 7 ) = - In Z EC + (Aec ~ A E c)o- (179) 

The validity of the mean-field approximation approach bases on the inequality (|178p that gives 
a low bound of the approximate free energy J^^Mh, g, 7). We determine the mean-field value 
Mt(g, r f) of the local mean-field action (|170p . which minimizes the approximate free energy (|179p 
and thus optimizes the low bound in Eq. (|178p . by satisfying the condition 

5 



SM b 



Jf c p (M h)5)7 ) 



= 0. (180) 

M h =M* 
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Using the mean-field value M£(g,j) and corresponding minimum of the approximate free energy 
J^$[M£(g, 7), g, 7] (|179p . we can gain some insights into the value of the 2-simplex area (|166|167p . 
and the critical points of the second-order phase transition, in terms of the gauge coupling g and 
Immirzi parameter 7. In addition, we can use the mean-field action (|170p with the value to 
calculate mean-field vacuum expectational values (• • -) to approximate true vacuum expectational 
values (• • •) that we discussed in Sees. IIIIHl IIII II and FVl 



B. Analytical calculations 

We can analytically calculate the mean- field partition function (|175p . First we integrate over 
quantized tetrad e M (x) and e u (x) fields, which is quadratic in Eq. (|170p (see Appendix [E]). Using 
the formula (|E2p . we have 

JJ f de^exp-AEC = Y[ det- l [I-T h ] (181) 

h€M heM 



and the Cayley-Hamilton formula for a determinant 



det _1 [J -T h ] = exp[-trln(/-r ?t )] 

= 1 + ^aa + 9 12^aaFbb + ^ab^ba) + ' 



2 

a a,b 

+ —\ 12 12 T aia Pl ■■■ T a n a Pn ( 182 ) 
• ai—a n p 

where P indicates permutations of (1, • • •, n) and Eq. (|182p is a sum of traces of symmetrized tensor 
products. The expression (|182p stops at the n-th order for a finite n-dimensional matrix T h in the 
space of the gauge group. 

Second we integrate over group- valued spin-connection U p (x + a p ), U p (x) and U v {x) fields 
defined at edges (x + a^,p), (x, fi) and (x, u) of the 2-simplex h(x) by using the properties of the 
invariant Haar measure: 

dU p {x) = 1 (183) 
dU^U^x) = (184) 
dU p {x)U; b (x)Ut cd {x') = ±-5^5 ac 6 bd 6(x-x'), (185) 

where dj = nj L nj R (rij Lt j R = 1jl,R + = 1/2,3/2,- • •) is the dimensions of irreducible 

representations j = (Jl, jr) of the gauge group SUl(2) <8> SUr(2), jn = ji = 1/2 and dj = 4 for 
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the fundamental representation. In Appendix [Ej we give more detailed calculations to obtain the 
mean-field partition function (|175|) . 

7 2 + 1 



Zec 



n 

h&M L 



1 + 



64 5 4 7 2 df 



(186) 



where Il/iex is the product of all 2-simplices h of the simplicial complex A4. The mean-field 
entropy is given by 

7 2 + l 



S = In Zec 



E ln 

h£M 



1 + 



A/" In 



where M = J2heM 1S ^ ne total number of 2-simplexes, and the mean- field free energy 



1 



\nZ EC 



-AHn 



1+ 6477^3 Mh ' 



(187) 



(188) 



where the inverse "temperature" /3 = 1/g 2 , see Eqs. (fT24l) and (fT29l) . 

We turn to calculate (Aec)o in Eq. (|178p . The mean-field value of Aec (|174p is calculated in 
Appendix [E] [see Eq. (|E7j) ]. 



= A/" ? 2 +A,M v 4 



1+ 7 \ + }^ m£ 



(189) 



32 9 4 7 2 (if 

where the vacuum expectational value with respect to the local mean-field partition function Zh 
(fT72l is defined by 

1 



<•••>? 



VUVe(- ■ -)e 



(190) 



The mean-field value (-4£c*)o (|189p has discrete values depending on the discrete values dj = 4, • • • 
of the fundamental state jl^r = 1/2 and excitation states jl,r = 3/2,- • •, coupling to different 
fermion spinor states ip L r- 

We are in the position to calculate (Aec)o in Eq. (|178p . Since there are three vertex fields 
in the smallest holonomy field X^v^U) (|120p that constitutes the regularized EC action Aec 
(|124p . (|129p . and (|13ip . while there is only one vertex field v VfJl in the mean- field action (|168p - 
(|170p . we assign the vertex field to the local mean-field action (|168H171]) of the 2-simplex h, the 
vertex-fields v pp ,v pu to the local mean-field actions of neighboring 2-simplices, and approximate 



(tr [v up {x)U p (x)v pp (x + a p )U p (x + a p )v pu {x + a u )U u ( 



x + a„)]) + h.c. 
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tr 



{vvpUnUpUyVupVpJ)* + h.c. 
« (^) 2 tr[(^C/ M C/ p ^)^ (t; w )J (<v)o] +h.c. 
^ (Zhftv^iv^UpU^ + h.c.) (v^) h {v pv ) h 

where (vp P v pi/ )^ = (v pu Vpp). Using Eqs. (|170p and (|17ip . we have 



{(tr 



>otr[<[eJ>* (M)^]}. 



(191) 



(192) 



In the last part of Appendix [Ej, we obtain 

1 ( 1 \ / 1 



(A. 



EC/o 



7 2 



Putting Eqs. (|187p . (| 189j) and f)193|) into the approximate free energy (|179p . we obtain 

.5/2 



ln(l + y) 



2y , y 



l + y 



i + y) 



where 



7 2 + l 



(193) 



(194) 



(195) 



In Fig. El we plot the approximate free energy (I179P as a function of the mean-field value Mh (|166p 
for selected values of the parameter \ (|195p . The minimal values of the approximate free energy 
3~EC (HZSJ) locate at the nonvanishing mean-field value M£ ^ 0, which increases as the parameter 
X decreases, namely, the gauge coupling increases. The gauge coupling g and Immirzi parameter 
7 remain to be determined. These two parameters (g, 7 ) should be determined at critical points of 
the second-order phase transition, as discussed in the last section. The mean-field approximation 
approach adopted here needs to be improved to see whether we can have a critical value Xc and for 
X > Xc the minimal value of the approximate free energy locates at the vanishing mean-field 
value ^ 0. It is usually difficult to study the vicinity of critical points of the second-order 
phase transition by the mean-field approximation approach. 

Considering the case that 7 > 1, dj = 4, g -» 4/3 for G cS ->■ G [see Eq. (|128p in Sec. IIIIF] , 
and x ~ 0.02, we have 



Mu* > 1, 



(196) 



see the curve for x = 0.03 in Fig. [UJ since becomes larger as x decreases. For larger gauge 
coupling g and higher dimensions dj of irreducible representations, the values of x (|195|) become 
smaller, and becomes larger. 
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FIG. 6: In the Planck unit a = 1, the approximate free energy (|179|) as a function of the mean-field value 
(|166[) is plotted for selected values \ — 0.03, 0.3, 3. The minimal values of the approximate free energy J 7 ^? 
locate at the nonvanishing mean- field value M£. The minimal locations are M£(x = 0.03) « 7.9, M£(x = 

0. 3)«2.1,MJ(x = 3)«0.8. 

Therefore, the mean-field value of the 2-simplex area (I166P 

(S h ) = a 2 M h * > a 2 = — ^— , (197) 

m Planck 

and the mean-field value of the volume element (|167j) 

<cfV(z)) = a 4 iV h (M h *) 2 > N h ^^. (198) 

m Planck 

Equations (|197|) and (|198p indicate that the averaged sizes of 2-simplex, 3-simplex, and 4-simplex, 

1. e., elements of the simplicial complex, are larger than the Planck length, which is probed by 
short wavelengths of quantum fields e^, Up, if) in strong gauge couplings g. This implies that due 
to the quantum gravity, the Planck length sets the scale for the minimal separation between two 
space-time points 25(. We end this section by noting that the mean-field approximation is not 



only a poor approximation, but also breaks diffeomorphism and local gauge symmetries. 
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VII. SOME REMARKS 



In addition to the Planck length a, the regularized EC action (|147j) proposed in this article 
contains three dimensionless parameters: the gauge coupling g; the Immirzi parameter 7 and the 
cosmological parameter A. In the view of the naive continuum limit, the regularized EC action (|147p 
proposed in this article is not unique. In principle, permitted by the diffeomorphism and local gauge 
invariances, the regularized action (I147j) is allowed to contain nonlocal high-dimensional (d > 6) 
operators of fields e^, Un and ip with extra free parameters. On the other hand, although the 
regularized EC action (|14Tj) approaches to the continuum EC action (|2ip in the naive continuous 
limit, it has not been clear yet whether the regularized EC theory is physically sensible. The 
regularized EC theory is physically sensible, only if only it has a nontrivial continuum limit, where 
we could possibly explore the relationship to the Minkowski counterpart. Therefore it is crucial, 
on the basis of nonperturbative methods and renormalization-group invariance, to find: 

1. the scaling invariant region (nontrivial ultraviolet fix points) (g c , j c , A c ), where the singu- 
larity in the free energy appears for phase transition occurring, and the physical correlation 
length £ of two-point Green-functions of fields is much larger than the Planck length, while 
the inverse correlation length gives the mass scale of low-energy excitations of the "ef- 
fective continuum theory" ; 

2. /3 function /3(g), i.e., the scale dependence of the gauge coupling g in the vicinity of the 
nontrivial ultraviolet fix points g c , and renormalization-group invariant equation 



in this scaling invariant region, and "constant" can only be obtained by nonperturbative 
methods. And it is a question how Eq. (I199p is related to 7 C and A c ; 

3. an effective action A e ^ c (|135p . all relevant and renormalizable operators [one-particle irre- 
ducible (1PI) functions] with effective dimension- four to obtain an effective low-energy theory 
in this scaling invariant region. 

The gauge-invariant correlation length £ can be possibly measured by the gauge-invariant two- 
point correlation function of the holonomy fields X^v, U) (|120p . 




(199) 



(X h [v(x),U(x)},Xl[v(y),U(y)}) 



~ e 



x-y\/i 



x- y\ > £, 



(200) 
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where \x — y\ indicates the separation between two holonomy fields Xh(v,U). Actually, Eq. (|200j) 
is related to the invariant curvature correlation function [see Eq. (lB12j) ], 

Although we have added the bare cosmological term (1139P into the regularized action, 1PI 
functions A e ^ c (|135p effectively contain this dimensional operator (|139p . which is related to the two- 
point correlation function (|200p . It is then a question what is the scaling property of this operator 
in terms of the low-energy scale £ -2 . We speculate that the gauge-invariant correlation length 
£, instead of the Planck length, sets the scale for the nonperturbative renormalized cosmological 
constant, i.e., 

Acosm ~ r 2 , (201) 

which is rather similar to the scale Aqcd calculated in the lattice QCD theory. This would possibly 
explain why the observed cosmological constant is much smaller than that expected in terms of the 
Planck scale [see Eq. (|199p ]. We also speculate that in the pure gravity at strong gauge coupling 
5 > 1, the scale £ -2 should measure the exponential area-decay law of holonomy fields (1134|138p 
for sufficiently large loops 

(X c (v, U)) ~ e' A — (c ^ 2 , A min (C) » f, (202) 
where A m [ n (C) is the minimal area, corresponding to the minimal number of 2-simplices h, that 



can be spanned by the loop C (see Ref. 26]). The scaling invariant region g c , scaling law (I199P and 
correlation length £ are important to study our present Universe (see Ref. [27I]). 

The effective quadralinear-fermion interactions in the continuum EC theory (|38|) are originated 
by integrating over static torsion fields and the torsion-free condition is satisfied as required by the 
equivalence principle. In this sense, quadralinear-fermion interactions are inevitable as long as the 
interacting between fermion and gravitational fields is included. 

The bilinear fermion action (|144p introduces a nonvanishing torsion field (I145D in the regularized 
EC theory. The torsion fields (I145P are not exactly static, however, they are fields only surviving 
in short distances at the Planck scale, which is due to the quantum gravity [see for example the 
mean-field approximation result (119611198]) ] . The effective quadralinear-fermion interactions (I146p is 
formulated by hand together with a torsion-free bilinear fermion action (I144p so that they approach 
to the fermion action of the continuum EC theory in the continuum limit. In principle, it should 
be possible to obtain an effective action by solving the discretized Cartan structure Eq. (|46p 
or Eq. (162p with the nonvanishing discretized torsion (I145p . and integrating over torsion fields 
at short distances, in the same way as (|30H38p of the continuum EC theory. In this way, one will 
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obtain a complicate effective action of fermion fields with high-order dimensional {d > 6) operators. 
However, we expect that in the continuum limit the relevant operators of fermion fields should be 
Eq. (I146P and its continuum counterpart (|40|) . 



On the other hand, due to the no-go theorem 



the bilinear fermion action (|144p has the prob- 



lem of either fermion doubling or chiral (parity) gauge symmetry breaking, which is inconsistent 
with the low-energy standard model for particle physics. As discussed, the effective quadralinear- 
fermion interactions (|146p are inevitable, due to mediating very massive torsion fields in short 
distances at the Planck scale. We expect that in the invariant scaling region of the nontrivial ul- 
traviolet fix points (g c , j c , A c ), the quadralinear-fermion interactions should be relevant operators, 



which not only give a possible resolution to the fermion doubling problem 
compelling dynamics for fermion mass generation 



33]. 
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30( | . but also the 



32| , via the Nambu Jona-Lasinio mechanism 
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Appendix A 

By using Eqs. JSgJ and ([MI) and the identity e A e^ = e A+B+[A,B}/2 ^ we ca i cu i ate JJ^x) (|57|) - 
in the native continuum limit: agu^ <C 1. Expanding U fJiU (x) in powers of agcu^, we have 

U^ v (x) = U^(x)U v (x + dp) = 

= exp liga[u^{x) +u v (x)] + iga 2 d^uj u (x) - ^(ga) 2 [u fl (x),uj l/ (x)} + 0(a 3 ) 



exp \ iga[uj^(x) +u v (x)] + iga 2 d^u v (x) - ~(ga) 2 [u ae (x) A w^a;)]^ + 0(a 3 ^ 



44 



= exp [igaa AB G AB + 0(a 3 )} , (Al) 

where 



= [^(x)+oj^(x)] + ad^(x) 
1 
2 



- ~(ga)lu M (x)Auj*(x)U, (A2) 



and 0(a 3 ) indicates high-order powers of agui^. In Eq. (jAl|) . we use [<J a b, a bc] = i$bb&ca (no sum 
with index 6) , [75 , a ca ] = and 

w Mi ,(x) = [upix^uivix)) = [u ae (x) Auj eb (x)}^ [a ae ,a eb ] = i[uj ae (x) A uj b e (x)} fW a ab . (A3) 

For exchanging //•<->■ v in Eqs. (jAip and ()A2p 



GiJf = [< B (*)+^(x)] + ad,< s (x) 



- iw[/(x)A We B (x)]„, (A4) 



As a result, the curvature R AB {x) ([19 



= a[d fl u AB {x)-d v u AB {x)] 



(ga)[u; A %x)Au B (x)U, (A5) 



where we use 



W Ae (x)A W f(x)]„ = -[uA(x) A W f (*)]„„. (A6) 



Appendix B 



The properties of the vertex fields v^ u (x) (|116p and (|117p : 



i 1 
V = 75^[7a7& -7&7a] 2 -(e£e i/ - e " e ^) 



2 

-750"afe(e a A e fe )^ = -v^ u = v Ufl (B2) 



for the case v /xu (x) = 75^ (x). Equations (jBip and ()B2p are the same for the case v^ u (x) = e Atl/ (x), 
because of 75 =75. For the sake of simplifying notations in following calculations, we introduce 
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+ab — -tab -tab — +ba arir \ „ — - , +ab 

L fiu b v\ii L fiv b iiv alLU z-^iv u ab L pu- 

We calculate the naive continuum limit of Eqs. (|120p . f)122[) . and (I123D in powers of gaujp. 
First, at the order 0(a ), we consider all link fields in Eqs. (|120p and (I122p to be identity, e.g., 
Up{x) ps l,U p (x + a^) ps 1, and U u (x + a„) ps 1. Using Eqs. (fT2!l - (fT23]l . (iBTTl . and (|B2|) . we obtain 
up to order O(a ) 

X h (v, U) + Xl(v, U) = tr [v V p{x)vp P {x + a„)v pv {x + a v )] + h.c. = 0. (B4) 

Second, at the order 0(a), we consider two link fields in Eqs. (|120p and (|122p to be identity. The 
case (1): U v (x + a v ) ps 1 and C/ p (x + a p ) ps 1, we have up to order 0(a), 

X h (v,U) ps tr [f^(x)[/ At (x)i; w (x + a^)v pv {x + a„)] 
ps tr K (U (ic)f (U p(2; + a fl )v pv (x + a„)] 

+ i 5 a^ B (x)tr[75CT a6 CTAija c df7 e/ ]^(x)t^ ) (x + a p )^(x + a^). (B5) 

for the case Vp V (x) = 75e At!/ (x). Using Eqs. (I120p - (I123I) . and (jB4j) . we have 

X ft («,E7) + 4(t;,E0 « i 5 a[^ B (a;)-^ B (x)]. 

fK(x)C(x + ap)t%(x + a u ). (B6) 

The case (2): Un[x + an) ps 1 and U p (x + a„) ps 1, we obtain the result with the replacement 
[u^ B (x) — uj„ b (x)\ — > [lj^ b (x) — u)^ B (x)] in Eq. (|B6|) . Taking into account all contributions from 
these cases, we obtain up to the order 0(a) 

X h (v,U)+xl(v,U) = 0. (B7) 

These results are also valid for the case Vpp(x) — 6 pi y(x), since the calculations of Eqs. ()B4p ~( jB6P 
without 75 are the same. 

Third, at the order 0(a 2 ), we consider one link field in Eqs. ()120p and (|122p to be identity, e.g., 
U p (x + ap) ps 1, 

X h (v,U) ps tT{v V p(x)Up(x)vpp(x + ap)v pv (x + a u )U u (x + a u )] 

ps tr[v V p(x)Up(x)U u (x)vpp(x + ap)vp U (x + a u )}, (B8) 

where in the second line, we use Eq. (|5o]) . [0^,75] = 0, [Up(x),v pu ] = 0(a), and U u (x + a v ) = 
U v (x) + 0(a). Using Eq. fl89|) or flS]) for f/^x) = Up(x)U v (x) and the result (jB"4) . we have up 
to 0(a 2 ) 

X h (v,U) ps tr [u I/A i(a;)J7 AtI/ (a;)u |t ^(a; + a^)v pv (x + a„)] 

= ia 5 G^ B (x)tr[7 5 CT a6 aAijo- c dCT e/ ]^(x)^' P (3; + a At )^(x + a^), (B9) 
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for the case v )11/ (x) = 75<7 Ati ,(x). Using the relationships X^iv, U) = Xh(y, £7)U<->k (|12ip and (jl22H 
and tf v = -tf p ([B3]), we have 

xl(v,U) « -iagG^{x)tx[^a ah a AB a cd a e} ]tll{x)tf p {x + a^t^x + a v ). (BIO) 

As a result, using Eq. (|A5P in Appendix^ we obtain up to 0(a 2 ) 

X h {v,U)+xl(v,U) « ia 2 5 ^(x)tr[7 5 a afc a AB a cd a e/ ]t^(x)t^(x + a M )^(x + a^). (Bll) 

For the case v pu (x) = e pu (x), the result is given by Eq. (|Blip without 75. 

In Appendix O we show the calculations of ^[^(TabCJAB^cd^ef] and tr[a a baABO'cd^ef] m 
Eq. (IBTT1) . Using these results (|C4]) and (lC8lh we obtain for the case iy«/(x) — 75^/w 

(x), 

X h (v, U) + X\(v, U) « 8a 2 5j R^(x)e aMB 4(x)t^(x + a„)t#(s + a,); (B12) 

and for the case v fMU (x) = e^ v (x), 

X h (v, U) + Xl(v, U) « 2i ■ 8a 2 gR^(x)t^(x)t^ p (x + a^tf^x + a v ). (B13) 

Using Eqs. (|119|) and (|B3p . we rewrite the fundamental area (|96|) and (j9?|) of the 2-simplex h(x) 
in terms of i2p(x + a^) and i*^(x + a v ): 

SU X + °m) = <r*S%(x + a M ), S£j(s + a M ) = -ia 2 ^ p (x + a„), (B14) 
^(x + c) = (7 cd 5^(x + a,,), S^(x + a„) = -ia 2 t cd {x + a*), (B15) 

where 5^ p (x + a M ) = — ^^(x + a^) and ^^(x+a^) = — S^x + a^). As discussed in Eqs. (|103p . ([96]) . 
and (f9T|) [see Sec. IIIIDj . the area of three area operators S pi/ (x), S p Jx + a M ) and S^ p (x + a v ) are 
identical. Therefore, equivalently to Eqs. (I104p and (1107p . we write the volume element contributed 
from the 2-simplex h{x) as 

dVh = Sf p {x + aix )S${x + a v ) = a% d p (x + a^)t c p d u (x + a u ) 

= S*(z)S$(z) = aH%(x)t%(x), (B16) 



where indexes c, d are summed, while indexes ^ and p are not summed. Using Eq. (|C7p in 
Appendix [Cj we obtain 

dU,(x) = S 2 (x) = 1 tr [^(x)5^( x )] , (B17) 

where 5^(x) = cr^S^x) and S#(x) = -ia?tf u {x). Using Eqs. ()B12IB13p and (IB14IIB17p . we 
can show the regularized Palatini action (|124|) and Host action (|129p approach to their continuum 
counterparts (f22|) and ([23]) in the naive continuum limit agu^ <C 1. 
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Appendix C 

It can be shown that tr[75<7 a b<7 c d(x e j] = for 75 = 70717273 hi the four-dimensional space-time, 
nonvanishing contributions of the following trace 

to[lS a abVABOcdOef\, (CI) 
come from the product of two spinor matrices cr's in Eq. (|C1|) being identical, 

tr [75 CTabVABVcdCref] => ^[l5^ab^AB] ■ (C2) 

In Eq. (jBllj) . as example, we take (i) a cc icr e f = 1 for c = e, d = f and (ii) a c d&ef = —1 c = /, d = e, 
E Wcd<ref]tfh(x + a^fplix + a u ) 

cdef 

= Y.^ a ^ cr cd]t c l f p (x + a^tf^x + a v ) + Y,[ a od^dc]tf p (x + a^)t%{x + a v ), 

cd cd 

= E O x + a ^( x + «-) - E O x + (* + a ») 

cd cd 

= 2 y £l%(x + aJ1%(x + a v ). (C3) 

cd 

Thus, in Eq. (jBllj) we have 



= 2tr[ 75 a a6 aA J B]C( a:: )^( :c + a M( x + a ») 



where we use the formula 



8ie abAB tl1(x)C(x + a,)C(x + a u ), (C4) 



tr 



{l^a AB ) = \ir ( 75 {a ab ,* AB }) = -4ie« bAB , (C5) 

and Eq. (|18|) . In the same way we calculate Eq. (jBllj) for other possibilities, e.g., cr a bO'ef = 1 f° r 
(i) a = e,b = f and (ii) (Tab&ef = —1 a = f,b = e. As a result, we obtain Eq. ()B12j) . 

Analogous to the discussions for Eq. (jC2j) . nonvanishing contributions to t^[(Tab^AB^cd^ef\ come 
from the product of two spinor matrices it's being identical, 

tr[a ab a A BcrcdVef} =^ tr[a ab a A B}- (C6) 

In Eq. (jBlip without 75, as example, we take (i) o c ^o e ^ = 1 for c = e, d = f and (ii) a c d<Jef = — 1 
c = f,d = e, and use formula 

trjVV 45 ) = A(5 aA 5 bB -5 aB 5 bA ). (C7) 
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As a result we obtain 

= 2tv[a ab a A B]t%{x)tf p (x + a„)tf v (x + a v ) 

= 2 • 8i£f (z)t°J(z + ap )t%(x + a u ), (C8) 
and Eq. ()Bli|) without 75 becomes Eq. ()B13j) , 



Appendix D 



Using the properties (|B1|) of the vertex field v fiv {x) = 75e Atl/ (x), we have 



X h (v,U) 



x\{v,U) 



e u (x)U p (x)U p (x + cip)U u (x + a u )ep(x) 



ep(x)Up(x)U p (x + cip)U u (x + a u )e u (x) 



Mi 



eJx)U v (x)U p (x + a v )Up(x + aje^a;) 



e u (x)U u (x)Up(x + a u )Up(x + a M )e A1 (x) 



M, 



hi 



(Dl) 



and 



Z ft ( U ,C/)+X^( U ,C7) 

= ^Mltrj 5 e v (x) Up(x)U p (x + ap)U u {x + a v ) - U v {x)U p (x + a u )Up(x + dp) e M (x) 

+^Mhtr7 5 e At (x) U v {x)U p {x + a v )U^{x + a^) - U p ,{x)U p (x + a^)U v {x + a v ) e u (x) 
= tr [e I/ (x)7 5 i7 I/M (x)e At (x)] - tr [eJx)^ 5 H u Jx)e u (x)] , 



(D2) 



where 75e^(x) = — e^(x)7s and the tensor 

H m (x) = l -Ml[U u {x)Up{x + a u )Up{x + ap)-Up{x)Up{x + ap)U u {x + a u )} 



-Mi 



U u (x)U p (x + a v )Up(x + dp) - UUx + cip)UUx + a v )Ul(x) 



-Ml [U u {x)U p {x + a u )Up{x + ap)} + h.c. 



-Mi 
2 h 



C/ i ,(x)?7 p (x + a I/ )^(x) 



+ h.a, 



(D3) 



H V p = -Hp V and H\ p = H U p, following the relations Up(x) = Up[(x + a M ), Ul(x) = U u (x + a u ) 
and U p (x + a^) = U^x + aj,). The Hp U (x) is a product of three edge fields U u (x), Up\(x) and 
U p (x + dp) of the 2-simplex h(x). For the case Vp U (x) = ep U (x), the same result can be obtained 
by the replacement 75 — > — 1 in Eq. (JD2]). The sum of two contributions gives Eqs. (I170p and (1171|) 
in the main text. 
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Appendix E 



For each 2-simplex h (fi ^ u ^ p), we have the fundamental area operator e M A e u = e v e^ — e^e v 
[see Eq. (|103p ] and tr(e v e M — e^e u ) = 0, we can rewrite the mean-field action fllTOj) as follows: 

A h = tr(e„e M - e^e v ) + A h 



tr 



(/-!%) 



tr < ( e„ e M ) 



where I is the identity matrix. For each single 2-simplex h, we have the integrations 



de^deu exp —Ah = det 1 [7 — T 71 ], 
de^e^exp-A = ^[-f — T^] "J det -1 [J — r' 1 ], 
1 de^e, e M „ exp = - T%1 -[I- r^det^/ - r* 



(El) 

(E2) 
(E3) 
(E4) 



Using Eqs. (|18ip and (|182p . we calculate the mean- field partition function (|175p 

= II / <ft/ M dc/ i ,d^det- 1 [i-r fe ] 



heM ' 



1 



n / du,du u du p i + £ v h aa + - Y,V h oA + rKj + 



(E5) 

heM - " 1 'a " a ,b 

In Eq. ()E5p . the first term is one due to the formula (|183p . the second term vanishes due to the 
formula (|184p . and nonvanishing contribution, due to Eqs. (|184p and (|185p . comes from the term 
r£ 6 r£ a in the third term. Using Eqs. (fTTTT) . <HM\ . and (fT85|) . we have 



1 



dU^dUp-^T^T 



h ■p'i 
ba 



1 ( 1 



a,6 



1/1,1 ' \) (t) j h dU ^ dU - dU p 



(75 - ^) .[^Wn[C/Jk(75 " [C^U[^]«[C4] 



1 ;' \ ) k 

2 V 8^2 / n 



1\ 2 



1 



1 



8? J + 7 



1 



tr 



(E6) 



As a result, we obtain the mean- field partition function (|186p in the main text. 

Using Eq. (IE4j) , we calculate the mean- field value of the mean-field action Ah (j!70j) of the single 
2-simplex fa, 



(A)o = (tr 



e T h e - e V h e 
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VUtv{ "" 



T h 



L j h wtr{2r^r^ + • • • jdet" 1 !/ - r h ] 



2Z h 

Z h \8g 2 
1/1 



Mc - -^tr 



1\2 
75 

7 



z h V8<? 2 

which gives Eq. (|189p in the main text. 
Using Eqs. ([E2]h (jE3jl . and (|E5|) and (T h 

i 1 



+ 1 



(E7) 



'MP 



-(T h ) pfl [see Eqs. (fT7il and ([D3])]. we have 



(W>J = ^/^^{[/-r^-fi-r^ldet- 1 !/-^ 

" 2r ^ + ...jdet- 1 [/-r h ] 



4 4 A 
i 2 

44 



1 



1\2 
75 

7 



(E8) 



and ([e pA ])5 = 



-([e,J)^. As a result, Eq. (H92J) becomes 



E 



(Zh 



{<tr 



e pFvp e v 



>otr 



7 2 + 1\ r/7 2 + 1\2 



(^) <E9) 



and we obtain Eq. (|193p in the main text. 
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